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ABSTRACT 

We present an overview of the Space Telescope A901/2 Galaxy Evolution Survey (STAGES). 
STAGES is a multiwavelength project designed to probe physical drivers of galaxy evolu- 
tion across a wide range of environments and luminosity. A complex multi-cluster system at 
z ~ 0.165 has been the subject of an 80-orbit F606W HST/ACS mosaic covering the full 
0.5° x 0.5° (^5x5 Mpc 2 ) span of the supercluster. Extensive multiwavelength observations 
with XMM-Newton, GALEX, Spitzer, 2dF, GMRT, and the 17-band COMBO- 17 photomet- 
ric redshift survey complement the HST imaging. Our survey goals include simultaneously 
linking galaxy morphology with other observables such as age, star-formation rate, nuclear 
activity, and stellar mass. In addition, with the multiwavelength dataset and new high resolu- 
tion mass maps from gravitational lensing, we are able to disentangle the large-scale structure 
of the system. By examining all aspects of environment we will be able to evaluate the relative 
importance of the dark matter halos, the local galaxy density, and the hot X-ray gas in driving 
galaxy transformation. This paper describes the HST imaging, data reduction, and creation 
of a master catalogue. We perform Sersic fitting on the HST images and conduct associated 
simulations to quantify completeness. In addition, we present the COMBO-17 photometric 
redshift catalogue and estimates of stellar masses and star-formation rates for this field. We 
define galaxy and cluster sample selection criteria which will be the basis for forthcoming 
science analyses, and present a compilation of notable objects in the field. Finally, we de- 
scribe the further multiwavelength observations and announce public access to the data and 
catalogues. 
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1 SURVEY MOTIVATION 

1.1 A multiwavelength approach to galaxy evolution as a 
function of environment 

The precise role that environment plays in shaping galaxy 
evolution is a hotly debated topic. Trends to passive and/or 
more spheroidal populations in dense environm e nts are widely 
observed: galaxy morphology (Dressier! Il980l : [Dressier et alJ 
19971; iGoto et al] 120031; iTreu et al.l 120031). co lour dKodama et ail 
200ll;lBlanton et a l. 2005; Baldr v et al.ll2006h . star-formation rate 
dGomez et al.ll2003l: iLewis et al.ll2003) . and stellar age and AGN 
fraction dKauffmann et al.ll2004h all correlate with measurements 
of the local galaxy densi ty. Furthermore, these relations persist over 
a wide r ange of re dshift ( Smith et al. 2005 ; Cooper et alj2007T) and 
density dBalogh et alj2004) . 

Disentangling the relative importance of internal and exter- 
nal physical mechanisms responsible for these relations is chal- 
lenging. It is natural to expect that high density environments will 
preferentially host older stellar populations. H ierarchical models 
of galaxy formation (e.g. lDe Lucia et al]|2006t) suggest that galax- 
ies in the highest density peaks started forming stars and assem- 
bling mass earlier: in essence they have a head-start. Simultane- 
ously, galaxies forming in high-density environments will have 
more time to experience the external influence of their local en- 
vironment. Those processes will also act on infalling galaxies as 
they are continuously accreted into larger haloes. There are many 
plausible physical mechanisms by which a gal axy could be trans- 
formed by its environment: removal of the hot dLarson et al.ll 1980l) 
or cold dGunn & Gottl 19721) gas supply throu gh ram-pressu re strip- 
ping; tidal effects leading to halo truncation (Bek kJl999l) or trig- 
gered star formation through gas compression (Fujita 1998); inter- 
actions betwe en galaxies themselves via low-speed major mergers 
dBarneslll992h or frequen t impulsive encounters termed 'harrass- 
ment' dMoore et al.fl998h . 

Though some of the above mechanisms are largely cluster- 
specific (e.g. ram-pressure stripping requires interaction with a hot 
intracluster medium), it is also increasingly clear that low density 
environm ents such as galaxy groups are important s ites for galaxy 
evolution dBalogh et al.ll2004l : IZabludoff et al.lll996l) . Additionally, 
luminosity (or more directly, mass) is also critical in regulating 
how suscept i ble a galaxy is to external influences. For example, 
lHaines et alJ d2006l) find that in low density environments in the 
SDSS the fraction of passive galaxies is a strong function of lu- 
minosity. They find a complete absence of passive dwarf galaxies 
in the lowest density regions (i.e., while luminous passive galaxies 
can occur in all environments, low-luminosity passive galaxies can 
only occur in dense environments). 

Understanding the full degree of transformation is further 
complicated by the amount of dust-obscured star formation that 
may or may not be p r esent. Many studies in the radio and MIR 
dMiller & Owenll2003l ; ICoia et alj|2005l ; lGallazzi et alj|2008[) have 
shown that an optical census of star formation can underestimate 
the true rate. Cluster-cluster variations are strong, with induced 
star formation linked to dynamically-disturbed large-scale struc- 
ture dGeach et alj|2006h . Nor are changes in morphology necessar- 
ily equivalent to changes in star formation. There is no guarantee 
that external processes causing an increase or decrease in the star- 
formation rate act on the same timescale, to the same degree, or 
in the same regime as those responsible for structural changes. A 
full census of star-formation, AGN activity, and morphology there- 
fore requires a comprehensive view of galaxies, including multi- 
wavelength coverage and high resolution imaging. These are the 



aims of the STAGES project described in this paper, targetting the 
Abell 901(a,b)/902 multiple cluster system (hearafter A901/2) at 
2 - 0.165. 

In addition to the STAGES coverage of A901/2, there are sev- 
eral other multiwavelength projects taking a similar approach to 
targeting large-scale structures. While we will argue below that 
STAGES occupies a particular niche, the following is a (non- 
exhaustive) list of surveys of large-scale structure including sub- 
stantial HST imaging. All are complementary to STAGES by way 
of the redshift range or dynamical state probed. The COSMOS 
survey has examined the evolution of th e morphology-density re- 
lation to z — 1.2 (Capaketal. 20 q3), paying partic ular atten- 
tion to a large s tructure at z = 0.7 I Guzzo et al.ll2007r) . Relevant 
to this work, in ISmolcic et al] (2007) they identify a complex of 
small clusters at 2 ~ 0.2 via a wide-angle tail radio galaxy. At 
intermediate redshift, an extensive comparison project has been 
undertaken targeting the two contrasting clusters CL0024+17 and 
MS0451-03 at 2 ~ 0.5 to compare the low- and high-luminosit y 
X-ray cluster environment (Moran et al. 20071; Icieach et al.ll2006h . 
Locally, the Coma cluster has also been extensively used as a 
laboratory f or galaxy evolution (Poggianti et al. 2004; Garter et al.l 
I2002l l2008). There are many other examples of cluster-focused en- 
vironmental studies covering a range of redsh ifts, including the 
large sample of EDisCS clusters at 2 > 0.5 dWhite et al J 1 20051 : 
IPoggianti etafll2006l ; iDesai et alj|2007h: and the ACS GTO cluster 
program of 7 clusters at 2 ~ 1 1 Postman et al.ll20o"5l : IGoto et aU 



l2005l : iBlakeslee etai]|2006l : iHomeier et alj2005h 

We summarize the motivation for our survey design as fol- 
lows. In order to successfully penetrate the environmental pro- 
cesses at work in shaping galaxy evolution, several areas must be 
simultaneously addressed: a wide range of environments; a wide 
range in galaxy luminosity; and sensitivity to both obscured and 
unobscured star formation, stellar masses, AGN, and detailed mor- 
phologies. Furthermore, it is essential to use not just a single proxy 
for 'environment' but to understand directly the relative influences 
of the local galaxy density, the hot ICM and the dark matter on 
galaxy transformation. A further advantage is given by examining 
systems that are not simply massive clusters already in equilibrium. 
By including systems in the process of formation (when extensive 
mixing has not yet erased the memory of early timescales), the var- 
ious environmental proxies listed above might still be disentangled. 

Therefore, the goal of STAGES is to focus attention on a single 
large-scale structure to understand the detailed aspects of galaxy 
evolution as a function of environment. While no single study will 
provide a definitive answer to the question of environment and 
galaxy evolution, we argue that STAGES occupies a unique van- 
tage point in this field, to be complemented by other studies locally 
and at higher redshift. 



1.2 Galaxy evolution as a function of redshift: STAGES and 
GEMS 

In addition to science focused on the narrow redshift slice contain- 
ing the multiple cluster system, the multiwavelength data presented 
here provide a valuable resource for those wishing to study the evo- 
lution of the galaxy population since 2 = 1. With the advent of the 
HST and multiwavelength data for this field, it is possible to quan- 
tify better the sample variance and investigate rare subsamples us- 
ing the combination of the STAGES field together with the Galaxy 
Evolution and Morphologies (GEMS; Rix et al]|2004h coverage of 
the Extended Chandra Deep Field South (CDFS). In particular, the 
HST data were chosen to have the same passband for both GEMS 
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(F606W and 850LP) and STAGES (F606W only, to allow study 
at optimum S/N of the cluster subpopulation and to optimise the 
weak lensing analysis). While the choice of F606W means that the 
data probe above the 4000A break for z < 0.5 only, for a num- 
ber of purposes the data can also be used at higher redshift (al- 
though in those cases one needs to be particularly cognizant of the 
effects of bandpass shifting and surface brightness dimming; such 
effects can be understood and calibrated using the GEMS 850LP 
and GOODS 850LP data). Furthermore, the 24[i.m observations 
(ED are well-matched in depth with the first Cycle GTO observa- 
tions of the CDFS ; analyses of the C DFS and A901/2 field s have 
been presented bv lZheng etall J2007h and lBell et ail d2007t) . Sev- 
eral projects are already exploiting this combined dataset (see ij5] 
for details), and with the publicly-available data in the CDFS, these 
samples provide a valuable starting point for many investigations 
of galaxy evolution. 

1.3 The Abell 901(a,b)/902 supercluster: a laboratory for 
galaxy evolution 

The A901/2 system is an exceptional testing ground with which 
to address environmental influences on galaxy evolution. Consist- 
ing of three clusters and related groups at z ~ 0.165, all within 
0.5° x 0.5°, this region has been the target of extensive ground- 
and space-based observations. We have used the resulting dataset 
to build up a comprehensive view of each of the main components 
of the large-scale structure: the galaxies, the dark matter, and the 
hot X-ray gas. The moderate redshift is advantageous as it enables 
us to study a large number of galaxies, yet the structure is contained 
within a tractable field-of-view and probes a volume with more gas 
and more star formation in general than in the local universe. 

The A901/2 region, centred at (q,<5)j 2 ooo = (9 h 56 m 17!3, 
— 10°01'11"), was originally one of three fields targeted by the 
COMBO- 17 survey dWolf et alj |2003). It was specifically chosen 
as a known overdensity due to the multiple Abell clusters present. 
These included two clusters (A901a and A901b) with X-ray lumi- 
nosities sufficient to be in cluded in the X-ray Brightest Abell-type 
Cluster Survey (XBACS; libeling et alj[l99a) of the R OSAT All- 
Sky Su rvey, though pointed ROSAT HRI observations by Schindler 
(2000) subsequently revealed that the emission from A901a suf- 
fers from confusion with several point sources in its vicinity. The 
extended X-ray emission in the field is further resolved by our 
deep XMM-Newton imaging (see i]4.6l >. Additional structures at 
z ~ 0.165 in the field include A902 and a collection of galaxies 
referred to as the Southwest Group (SWG). 

The five broad- and 12 medium-band observations from 
COMBO- 17 provide high-quality photometric redshifts and spec- 
tral energy distributions (SEDs). Together with the high-quality 
imaging for ground-based gravitational lensing, the A901/2 data 
have been used in a variety of papers to date. COMBO-17 derived 
results include 2D and 3D reconstruc tions of the mass distribution 
dGrav et alj|2002l; Ijavlor et alj [20041) : the star-formation-density 
relation iGrav et al.l2004h : the discovery of a subs tantial population 
of intermediate-age, dusty red cluster galaxies (Wolfetal. 2005, 
here- after WGM05); and the mor phology-density dLane et ail 
2007) and morphology-age-density dWolf et al J 20 07?) relations. 

Further afield, the clusters are also known to be part of a larger 
structure together with neighbouring clusters Abell 907 and Abell 
868 (1.5 degrees and 2.6 degrees away, respectively). Nowak et al. 
(in prep.) used a percolation (also called ' friends-of-friends') al - 
gorithm on the REFLEX cluster catalogue (Bohringer et alj|2004h 
to produce a catalogue of 79 X-ray superclusters. Entry 33 is the 



A868/A901a/A901b/A902/A907 supercluster, which also contains 
an additional, but not very bright, non-Abell cluster. Though not 
observed as part of the STAGES study, these clusters are included 
in the constrained N-body simulations used to understand the for- 
mation history of the large-scale structure (i ]4.8l >. 

The plan of this paper is as follows: in Sj2] we outline the ob- 
servations taken to construct the 80-tile mosaic with the Advanced 
Camera for Surveys on HST We discuss data reduction, object de- 
tection, and Sersic profile fitting. In Sj3]we present the COMBO-17 
catalogue for the A901/2 field and discuss how the two catalogues 
are matched. In Sj4]we present a summary of the further multiwave- 
length data for the field and derived quantities such as stellar masses 
and star-formation rates. We finish with describing ongoing science 
goals, future prospects, and instructions for public access to the 
data and catalogues described within. Appendix [A] contains details 
on ten individual objects of particular interest within the field. 

Throughout this paper we adopt a concordance cosmology 
with Q. m = 0.3, n x = 0.7, and H = 70 km s" 1 Mpc" 1 . 
In this cosmology, 1" = 2.83 kpc at the redshift of the super- 
cluster (z ~ 0.165), and the COMBO-17 field-of-view covers 
5.3 x 5.1 Mpc 2 . Magnitudes derived from the HST imaging (ij2} 
in the F606W (V-band) filter are on the AB systemQ while magni- 
tudes from COMBO-17 (® in all filters are on the Vega system. 



2 HST DATA 
2.1 Observations 

The primary goal of the STAGES HST imaging was to obtain mor- 
phologies and structural parameters for all cluster galaxies down to 
R = 24 (M v ~ -16 at z ~ 0.165). The full area of the COMBO- 
17 observations was targeted to sample a wide range of environ- 
ments. Secondary goals included obtaining accurate shape mea- 
surements of faint background galaxies for the purposes of weak 
lensing, and measuring morphologies and structural parameters for 
all remaining foreground and background galaxies to R = 24. 
As discussed in 1 .31 the survey design and filte r was chosen to 
match that of the GEMS survey dRix et al.|[2004h of the Chandra 
Deep Field South (CDFS). The CDFS is another field with both 
COMBO-17 and HST coverage, but in contrast to the A901/2 field 
is known to contain little significant large-scale structure. It will 
therefore serve as a matched control sample for comparing cluster 
and field environments at similar epochs. 

To this end we constructed an 80-tile mosaic with ACS in Cy- 
cle 13 to cover an area of roughly 29.5' x29. 5' in the F606W filter, 
with a mean overlap of 100 pixels between tiles. Scheduling con- 
straints forced the roll angle to be 125 degrees for the majority of 
observations, and one gap in the northeast corner was imposed on 
the otherwise contiguous region due to a bright (V = 9) star. A 4- 
point parallelogram-shaped dithering pattern was employed, with 
shifts of 2.5 pixels in each direction. An additional shift of 60.5 pix- 
els in the y-direction was included between dithers two and three 
in order to bridge the chip gap. 

Concerns about a time- varying PSF and possible effects on the 
weak lensing measurements drove the requirement for the observa- 
tions to be taken in as short a time frame as possible. In practice 
this was largely successful, with > 50% of tiles observed in a sin- 
gle five-day period (Fig. [TJ, and > 90% within 21 days. Six tiles 
(29, 75, 76, 77, 79, 80) were unobservable in that cycle and were 



For F606W, mAB — m~Ve 
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Figure 1. Cumulative plot of ACS data acquisition. In order to minimize 
the effects of a time-vary PSF on weak lensing applications, 50% of tiles 
were taken within 5 days and 90% within 21 days. The remaining 7 tiles 
were observed 6 months later. 



re-observed six months later, with a 180 degree rotation. Further- 
more, tile 46 was also re-observed at this orientation as the original 
observation failed due to a lack of guide stars. These seven tiles 
were observed following the transition to two-gyro mode with no 
adverse consequences in image quality. 

Details of the observations are listed in Table Q] A schematic 
of the field showing the ACS tiles and the multiwavelength obser- 
vations is shown in Fig. [2] Additionally, four parallel observations 
with WFPC2 (F450W) and NICMOS3 (Fl 10W and F160W) were 
obtained simultaneously for each ACS pointing. Due to the separa- 
tion of different instruments on the HST focal plane, most but not 
all parallel images overlap with the ACS mosaic (52/10/18 WFPC 
images and 42/9/29 NICMOS3 images have full/partial/no overlap 
with the ACS mosaic; most NICMOS3 images have partial over- 
lap with a WFPC2 image). In this paper we restrict ourselves to a 
discussion of the primary ACS data, analysis of the parallels will 
follow in a future publication. 



2.2 ACS data reduction 

We retrieve the reduced STAGES images processed by the CAL- 
NICA pipeline of STScI, which corrects for bias subtraction and 
flat-fielding. However, as the ACS camera is located 6 arcmin off 
the centre of the HST optical axis, the images from the telescope 
have a field-of-view with a parallelogram keystone distortion. To 
produce a final science image from the reduced pipeline data, we 
therefore also have to remove the geometric distortion before com- 
bining the individual dithered sub-exposures. The removal of the 
image distortion is no w fairly routine through t he use of the MUL- 
TIDRIZZLE software ( Koekemoer et al . 2007). However, our par- 
ticular science goals motivated us to make several changes when 
optimizing the default settings and combining the raw images. 
These changes are discussed below. 



2.2.1 Image Distortion Correction 

In STAGES, the science driver that demands the highest quality 
data reduction in terms of producing the most consistent and sta- 
ble PSF from image to image, and across the field of view, is weak 



lensing (Hevmans et al. 20 081). With this goa l in mind, we bene- 
fit from the experience of lRhodes et al. (2007), who conducted de- 
tailed studies of how the pixel values are re-binned when the images 
are corrected for image distortion. Briefly speaking, to transform 
an image that is sampled on a geometrically distorted grid onto 
one that is a uniform Cartesian grid fundamentally involves rebin- 
ning, i.e. interpolating, the original pixel values into the new grid. 
Doing so is not a straightforward process since the original ACS 
pixel scale samples the telescope diffraction limit below Nyquist 
frequency, i.e. the telescope PSF is undersampled. When a PSF 
is undersampled, aliasing of the pixel fluxes occurs, the result of 
which is that the recorded structure of the PSF appears to change 
with position, depending on the exact sub-pixel centroid of the PSF. 
This variability effectively produces a change in the ellipticity of 
the PSF as a function of sub-pixel position, even if the PSF should 
be identical everywhere. Because stellar PSFs are randomly cen- 
tred about a pixel the intrinsic ellipticity one then measures has a 
non-zero scatter. So, as weak lensing relies heavily on measuring 
the ellipticities of galaxies, which are convolved by the PSF, the 
scatter in the PSF ellipticity contributes significant noise to weak 
lensing measurements. 

An additional issue with non-Nyquist sampled images is that 
the process of interpolating pixel values necessarily degrades the 
original image resolution. While the intrinsic resolution can in prin- 
ciple be recovered by dithering the images while making observa- 
tions, strictly speaking this inversion is only possible when the im- 
age is on a perfect Cartesian grid at the start, i.e. with no image dis- 
tortion. Otherwise, there would be a residual "beating frequency" 
in the sampling of the reconstituted image, such that some pixels 
would be better sampled than others. Because of this, recovering 
the intrinsic resolution of the telescope when the field is distorted 
is not a well posed problem, and cannot easily be solved by a small 
number of image dithers. Some resolution loss will necessarily oc- 
cur in some parts of the image. This is especially true if the final 
images are combined after having been geometrically corrected, as 
is currently the process in MULTIDRIZZLE. One last, unavoid- 
able, side effect of interpolating a non-Nyquist sampled image is 
that the pixel values become necessarily correlated. However, the 
degree of resolution loss and noise correlation can be balanced by a 
suitable choice of interpolation kernels: whereas square top-hat ker- 
nels effectively amounts to linear interpolation and correlates only 
the immediate neighbour pixels but cause high interpolation (pixel- 
lation) noise, bell-shaped kernels (e.g. Gaussian and Sine) correlate 
more pixels but better preserve the image resolution. 

In light of these issues, it is clear that the goal of an optimal 
HST data reduction should be a dataset where the PSF structure is 
stable across the field of view and reproducible from image tile to 
tile. The contribution to the PSF variation by the stochastic aliasing 
of the PSF that necessarily occurs during 'drizzling' can be reduced 
by appropriate choi ces of drizzling kernel and output pixel scale. 
Rhodes etalJ ( l2007h characterize PSF stability in terms of the scat- 
ter in the apparent ellipticity of the PSF in the ACS field of view. 
After experimenting, they determine that the optimal set of param- 
eters in MULTIDRIZZLE to use is a Gaussian drizzling kernel, 
pixf rac=0.8, and an output pixel scale of 0'.'03. We thus follow 
their approach by adopting those parameters for our own reduction, 
while keeping all the other default parameters unchanged. How- 
ever, they note, as we do, that a Gaussian kernel causes more corre- 
lated pixels than tophat kernels. Nonetheless because the choice of 
interpolation kernel amounts effectively to a smoothing kernel, cor- 
related noise should in principle not have an impact on photometry 
statistics since the flux is conserved. Moreover, the same interpola- 
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Date 


a 


6 


Exposure 


N h ot 


N co ld 
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[dd/mm/yyyy] 


[J2000] 


[J2000] 


[s] 








1 


09 07 2005 


09:55:22.8 


-10:14:01 


1960 


85 1 


173 


796 


2 


07 07 2005 


09:55:44.5 


-10:13:54 


1960 


1082 


209 


982 


3 


08 07 2005 


09:55:33.4 


-10:12:03 


1960 


1157 


233 


1008 


4 


07 07 2005 


09:55:22.4 


-10:10:06 


1960 


1051 


199 


927 


5 


04 07 2005 


09:56:09.5 


-10:14:26 


1950 


1069 


195 


973 


6 


03 07 2005 


09:55:58.7 


-10:12:33 


1950 


1151 


219 


1027 


7 


04 07 2005 


09:55:47.9 


-10:10:39 


1950 


1038 


237 


905 


8 


04 07 2005 


09:55:37.0 


-10:08:45 


1950 


1095 


262 


938 


9 


04 07 2005 


09:55:26.2 


-10:06:52 


1950 


1020 


188 


876 


10 


05 07 2005 


09:55:15.4 


-10:04:58 


1950 


1014 


184 


938 


1 1 


07 07 2005 


09:56:38.6 


-10:15:34 


1960 


989 


219 


876 


12 


04 07 2005 


09:56:27.8 


-10:13:40 


1950 


1020 


226 


885 


13 


28 06 2005 


09:56:16.9 


-10:11:46 


2120 


1193 


256 


1037 


14 


28 06 2005 


09:56:06.1 


-10:09:53 


2120 


1391 


254 


1111 


15 


28 06 2005 


09:55:55.3 


-10:07:59 


2120 


1182 


253 


1052 


16 


29 06 2005 


09:55:44.5 


-10:06:06 


2120 


1109 


208 


940 


17 


29 06 2005 


09:55:33.7 


-10:04:12 


1960 


1116 


250 


888 


18 


04 07 2005 


09:55:22.9 


-10:02:18 


1950 


995 


178 


868 


19 


09 07 2005 


09:56:57.3 


-10:14:25 


1960 


963 


180 


786 


20 


07 07 2005 


09:56:46.0 


-10:12:54 


1960 


979 


222 


829 


21 


30 06 2005 


09:56:35.2 


-10:11:00 


1960 


1166 


288 


1005 


22 


28 06 2005 


09:56:24.4 


-10:09:07 


2120 


1 193 


263 


1012 


23 


25 06 2005 


09:56:13.6 


-10:07:13 


2120 


1143 


241 


1000 


24 


25 06 2005 


09:56:02.8 


-10:05:19 


2120 


1244 


254 


1128 


25 


22 06 2005 


09:55:52.0 


-10:03:26 


2120 


1274 


248 


1051 


26 


29 06 2005 


09:55:41.1 


-10:01:32 


1960 


1214 


275 


1063 


27 


05 07 2005 


09:55:30.3 


-09:59:39 


1950 


1258 


279 


1068 


28 


08 07 2005 


09:55:19.5 


-09:57:45 


1960 


1161 


220 


1052 


29 


04 01 2006 


09:57:10.7 


-10:14:08 


2120 


1274 


272 


1123 


30 


09 07 2005 


09:57:04.5 


-10:11:48 


1960 


943 


209 


781 


31 


08 07 2005 


09:56:53.5 


-10:10:14 


1960 


900 


200 


713 


32 


03 07 2005 


09:56:42.7 


-10:08:20 


1950 


1023 


214 


884 


33 


28 06 2005 


09:56:31.9 


-10:06:27 


2120 


1150 


223 


955 


34 


22 06 2005 


09:56:21.0 


-10:04:33 


2120 


1318 


243 


1111 


35 


22 06 2005 


09:56:10.2 


-10:02:40 


2120 


1220 


244 


1028 


36 


24 06 2005 


09:55:59.4 


-10:00:46 


2120 


1320 


287 


1101 


37 


29 06 2005 


09:55:48.6 


-09:58:53 


1960 


1150 


239 


974 


38 


05 07 2005 


09:55:37.8 


-09:56:59 


1950 


1123 


205 


951 


39 


08 07 2005 


09:55:27.0 


-09:55:05 


1960 


1094 


210 


965 


40 


09 07 2005 


09:57:12.5 


-10:09:14 


1960 


1062 


198 


916 


41 


07 07 2005 


09:57:00.9 


-10:07:34 


1960 


962 


176 


828 


42 


03 07 2005 


09:56:50.1 


-10:05:41 


1950 


1090 


205 


928 


43 


27 06 2005 


09:56:39.3 


-10:03:47 


2120 


1198 


202 


1052 


44 


27 06 2005 


09:56:28.5 


-10:01:54 


2120 


1266 


230 


1046 


45 


23 06 2005 


09:56:17.7 


-10:00:00 


2120 


1280 


285 


1064 


46 


01 01 2006 


09:56:05.4 


-09:57:47 


2120 


1438 


355 


1235 


47 


01 07 2005 


09:55:56.0 


-09:56:13 


1960 


1198 


273 


972 


48 


06 07 2005 


09:55:45.2 


-09:54:19 


1950 


989 


176 


852 


49 


06 07 2005 


09:55:34.4 


-09:52:26 


1960 


1054 


223 


901 


50 


09 07 2005 


09:55:24.4 


-09:50:31 


1960 


984 


212 


832 


51 


07 07 2005 


09:57:08.4 


-10:04:55 


1960 


1050 


189 


923 


52 


03 07 2005 


09:56:57.6 


-10:03:01 


1960 


1142 


209 


941 


53 


03 07 2005 


09:56:46.8 


-10:01:07 


1950 


1135 


21 1 


920 


54 


02 07 2005 


09:56:36.0 


-09:59:14 


1950 


1131 


228 


921 


55 


02 07 2005 


09:56:25.1 


-09:57:20 


1960 


1205 


31 1 


974 


56 


02 07 2005 


09:56:14.3 


-09:55:27 


1960 


1097 


242 


891 


57 


01 07 2005 


09:56:03.5 


-09:53:33 


1960 


1090 


210 


911 


58 


06 07 2005 


09:55:52.7 


-09:51:40 


1950 


1130 


201 


975 


59 


08 07 2005 


09:55:32.7 


-09:48:15 


1960 


1075 


204 


900 


60 


07 07 2005 


09:57:15.8 


-10:02:15 


1950 


1028 


183 


912 



Table 1. Details of STAGES HST/ACS observations. Only the second (successful) acquisition of tile 46 is listed. 'Hot', 'cold', and 'good' SExtractor configu- 
rations are described in £|2.3I Tiles 29, 46, 75, 76, 77, 79, and 80 are oriented at 180° with respect to the rest of the mosaic. The exposure time varied according 
to the maximum window of visibility available in each orbit. 
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Table 1 - continued 



Tile 


Date 


a 


5 


Exposure 


Nhot 


N C old 


Ngood 




[dd/mm/yyyy] 


[J2000] 


[J2000] 


[s] 








01 


m m inns 
u / u / ZUUJ 


no-s7-ns n 
uy. j / .uj.u 


i o-nn-9 1 

-1U.UU. Zl 


1 QSO 

iyju 


07 1 


1 R 3 
1 OJ 


R9A 
oZO 


A9 

oz 


m (T7 inns 

U / U / ZUUJ 


no-s^-s/i i 

UV. JO. J4.Z 


no-ss-9fi 

-Uy. JO.ZO 


1 osn 
iyju 


1 ns9 

IUjZ 


1 R 1 
1 o4 


on i 


A3 

Oj 


A A 07 9 fins 

UO U / ZUUJ 


Uy .JO. 4J.4 


OQ-SA-^A 

-uy. jo. J4 


1 QSO 

iyju 


1 1 A1 
1141 


9 1 7 
Z 1 / 


y ju 


04 


HA (Y7 7fV")S 
UO U/ ZUUJ 


no-s^-^i a 
uy. jo. jz.o 


no-s/i -a i 
-uy. 04.41 


1 asn 
iyju 


1 n^o 
luoy 


9 99 

ZZZ 


Ron 


fSS 


t\f. 07 inns 

uo U / ZUUJ 


OQ-SA-91 R 

uy . jo.z i .o 


OQ-S9-A7 

-uy . jz.4 / 


1 QSO 

iyju 


1 07 1 
1U / 1 


997 

ZZ / 


GOR 

yuo 


00 


cxfi, m inns 

UO U/ ZUUJ 


no ■ s £, ■ 1 1 n 
uy. jo. i i .u 


no-sn-s^i 
-uy. ju. j j 


1 QSO 

iyju 


1 m a 

1U14 


999 

zzz 


RSO 

o jy 


o / 


n^ 07 inns 

UO U / ZUUJ 


OQ-s^-on 1 
uy . jo.uu. i 


-uy .4o.OU 


1 QSO 

iyju 


1 OAA 
1U40 


99A 
ZZO 


Q91 

yzz 


A.R 
Do 


nc m inns 

Uo U/ ZUUJ 


no-ss-/io ^ 
uy. j j.4y. j 


no-/i7-nA 
-uy.4 / .uo 


1 QAO 

lyou 


Q£7 

yo / 


1 7Q 

i /y 


RS 1 


ov 


1 07 inns 

1 U U / ZUUJ 


OQ-S7- 1 7 s 

uy. j / . iz. j 


OQ-S7-A9 

-uy . j / .4Z 


1 QAO 

lyou 


R7A 
o /O 


1 AS 
1 4 J 


7RA 
/ o4 


/U 


no m inns 
uy U/ ZUUJ 


no-s7-m 7 
uy . j / . u i . / 


no-ss-/i r 
-uy. jj.4o 


1 Q^O 

lyou 


O 1 -1 

yj4 


1 R ~K 
1 OJ 


70R 


1 1 


no 07 ions 
uy u / zuuj 


uy. jo. ju.y 


OQ-S^-SA 

-uy. j j. j4 


1 QAO 

lyou 


1 0^9 
IUJZ 


1 R1 
1 o_ 


( N . a 


79 

1 z 


1 o 07 ions 

1 u U / ZUUJ 


OQ-S^-AO n 

uy . jo.'tu.u 


OQ-S9-0 1 

-uy. jz.ui 


1 Q£0 

lyou 


1 1 1 R 
1 1 1 o 


1 1 1 


yju 


7^ 
/ J 


no 07 ions 

\)y U / ZUUJ 


OQ-SA-9Q 7 

uy . jo.zy.z 


0Q-S0-07 

-uy . ju.u / 


1 QAO 

iyou 


Q 1 O 

y i u 


1 /SR 
1 Do 


77"^ 
ill 


74 


08 07 2005 


09:56:18.4 


-09:48:14 


I960 


907 


192 


822 


75 


04 01 2006 


09:57:11.0 


-09:53:30 


2120 


1708 


260 


1140 


76 


05 01 2006 


09:57:00.3 


-09:51:39 


2120 


1444 


275 


1134 


77 


05 01 2006 


09:56:49.5 


-09:49:48 


2120 


1324 


287 


1094 


78 


05 07 2005 


09:56:40.6 


-09:48:11 


1960 


1031 


184 


842 


79 


05 01 2006 


09:57:12.9 


-09:50:05 


2120 


1357 


302 


1019 


80 


05 01 2006 


09:57:02.8 


-09:48:36 


2120 


1255 


246 


973 



tion (smoothing) kernel propagates into the PSF, thus the choice of 
kernel should also not impact galaxy fitting analyses. 

2.2.2 Sky pedestal and further image flattening correction 

The images obtained from the HST archive have been bias sub- 
tracted and flatfielded. However, large-scale non-flatness on the or- 
der of 2-4% remains in the images, and there are slight but notice- 
able pedestal offsets that remain between the four quadrants. These 
large scale patterns and pedestals are both stationary and consis- 
tent in images that are observed closely in time. And even though 
MULTIDRIZZLE tries to equalize the pedestals before combin- 
ing the final images, the correction is not always perfect due to 
object contamination when computing the sky pedestal. These ef- 
fects are small, and the sky pedestal issue only affect large objects 
situated right on image boundaries, so that the effects on the en- 
tire survey itself may only be cosmetic. Nevertheless, we try to 
correct for the effects by producing a median image of data ob- 
served closely in time, after first rejecting the brightest 30% and 
faintest 20% of the images (to avoid over-subtraction). Then, for 
each of the four CCD quadrants, we fit a low order 2-D cubic- 
spline surface (IRAF/imsurfit) individually to model the large scale 
non-uniformity in the median sky image, and to remove noise. The 
noiseless model of the sky is then subtracted from all the data ob- 
served closely in time. After correction, the mean background in the 
four quadrants is essentially equal, and the residual non-flatness is 
< 1%. 



2.3 Object Detection 

Object detection and cataloguing were carried out automatically 
on the S TAGES F606W imaging data using the SExtractor V2.5.0 
software dBertin & Arnouts 1996). An optimized, dual ('cold' and 
'hot') configuration was used, following the strateg y developed 
for H ST/ACS data of similar depth for GEMS dCaldwell et all 
2008). The main challenge to extracting sources from the STAGES 



ACS data is the tradeoff between deblending high-surface bright- 
ness cluster members that are close on the sky in projection, and 
avoiding spurious splitting ('shredding') of highly structured spi- 
ral galaxies into multiple sources. In addition, we desire high de- 
tection completeness for faint, and often low-surface brightness, 
background galaxies. To optimize the detection completeness and 
deblending reliability for counterparts to i? ap ^ 24 mag galaxiefl 
from the COMBO- 17 catalogue, we fine-tuned the combination of 
cold and hot configuration parameters using three representative 
STAGES tiles (21, 39, and 55). For STAGES, we converged on the 
parameters given in Table [2] which successfully detected 99.5% 
(650/653) of the 7? ap sC 24 mag COMBO- 17 galaxies on these 
tiles, with reliable deblending for 98.0%. 

SExtractor produces a list of source positions and basic pho- 
tometric parameters for each astrometrically/photometrically cali- 
brated image, and produces a segmentation map that parses the im- 
age into source and background pixe ls, which is neces sary for sub- 
sequent galaxy fitting with GALFIT jPeng et al.l2002h described in 
£)2.4| For both configurations, a weight map (cx variance -1 ) and a 
three-pixel (FWHM) top-hat filtering kernel were used. The former 
suppresses spurious detections on low-weight pixels, and the latter 
discriminates against noise peaks, which statistically have smaller 
extent than real sources as convolved by the instrumental PSF. Our 
final catalogue contains 75 805 unique F606W sources uniformly 
and automatically identified from 17 978 objects detected in the 
cold run, and 89 464 'good' sources found in the hot run (before 
rejection of the unwanted hot detections that fell within the isopho- 
tal area of any cold detection). A total of 5 921 objects were man- 
ually removed from the catalogue after the detection stage. These 
detections are mainly over-deblended galaxies or image defects like 
cosmic rays. Another set of 658 detections were included in fitting 
the sample galaxies to ensure the accurate fitting of real objects, 
but excluded from the final catalogue. These were also mainly cos- 



2 COMBO-17 redshifts are mostly useful at i? ap $C 24 for reasons dis- 
cussed in detail in [fj] and so we adopt this cut for our main science sample. 
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Figure 2. Layout of multiwavelength observations of the A901/2 field. The numbered tiles represent the 80-orbit STAGES mosaic with HST/ACS, which 
overlaps the 31.5x30 arcmin COMBO-17 field-of-view (long-dashed square). The seven shaded tiles were observed ~6 months after the bulk of the obser- 
vations and with a 180° rotation. The centres of A901a/A901b/A902/SWG are found in tiles 55/36/21/8 respectively. Interior to the STAGES region are the 
XMM-Newton coverage (heavy solid polygon) and the GMRT 1280 MHz observations (short-dashed circle, indicating half-power beam width). The STAGES 
area is also overlapped by the field-of-view of the Spitzer 24|am imaging (solid polygon), the GMRT 610 MHz observations (long-dashed circle), and the 
GALEX imaging (dotted circle). 



mic ray hits or stellar diffraction spikes. Although the main analysis 
was performed on a tile-by-tile basis, rather than mosaic-wise, the 
main catalogue only contains unique sources. Objects detected on 
two tiles enter the catalogue only once. The most interior-located 
was selected for entry into the catalogue. The breakdown of cold, 
hot, and good sources per ACS frame is given in Table[T] 

In Fig. [3] we show a histogram of various object samples in 
the region of the HST-mosaic that overlaps with COMBO-17. The 
HST data start to become incomplete at V606 ~ 26 (solid line). 
Stars (hashed histogram) only make up a significant fraction of all 
detections at the brightest magnitudes. A histogram of counterparts 
from a cross-correlation with COMBO-17 is shown in light grey. 
When the match is restricted to extended objects with i? ap < 24 
(ie. the primary 'galaxy' sample for which we have reliable photo- 
metric redshifts), the HST sources largely have Vt>06 < 24. 

Star-galaxy separation is performed in the apparent magnitude 
- size plane spanned by the SExtractor parameters MAG_BEST 
(Veoe) and FLUX_RADIUS (r/). Objects with 

log(r/) < max (0.35; 1.60 -0.05^606; 5.10 -0.22\4oe) (1) 
are classified as point sources; sources above that line are identified 



as extended sources (galaxies). This plane is shown in Fig. [4] The 
separation line clearly delineates compact and extended sources, in 
particular when inspecting the COMBO-17 sources only (crosses). 
Note that those AGN for which the point source dominates are also 
found on the point- source locus and therefore are removed from the 
galaxy sample by this selection. 

In the Fig. [5] we display the galaxy fraction as a function of 
Veoe magnitude (grey histogram). Out to Veos ~ 22 almost every 
galaxy detection on the HST images has a COMBO-17 counterpart; 
at the COMBO-17 sample limit V606 ~ 24 the matching complete- 
ness for STAGES objects is still ~90%. The cross-matching be- 
tween COMBO-17 and the HST data is described in more detail 
in § 13.21 where completeness is defined in reverse, i.e. maximizing 
HST counterparts for COMBO-17 objects. 



2.4 Sersic profile fitting 

To obtain Sersic model fits for each STAGES galaxy, the imaging 
data were processed with the data pipeline GALAPAGOS (Galaxy 
Analysis over Large Areas: Parameter Assessment by GALFITting 
Objects from SExtractor; Barden et al., in prep.). GALAPAGOS 
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Table 2. Dual SExtractor parameter values for STAGES F606W object detection in 'cold' and 'hot' configurations. 



Parameter 


Cold 


Hot 


Description 


DETECT.THRESH 


2.8 


1.5 


detection threshold above background 


DETECTJVIINAREA 


140 


45 


minimum connected pixels above threshold 


DEBLENDJVIINCONT 


0.02 


0.25 


minimum flux/peak contrast ratio 


DEBLENDJSfTHRESH 


64 


32 


number of deblending threshold steps 



7000 



6000 r 




16 18 20 22 24 26 28 
V 6ob [mag] 



Figure 3. Source detections in the HST mosaic (overlap region with 
STAGES and COMBO-17 coverage). The solid line represents all SExtrac- 
tor detected sources (74 534 objects). The grey histograms shows all objects 
with a corresponding match in the COMBO- 1 7 catalogue (light grey; 50 701 
sources) and extended sources with i?ap < 24 (dark grey; 12 748 sources). 
In addition, the hashed region indicates stars as defined by our star-galaxy 
separation criterion (Equation [T] 4969 stars in total). In the inset we high- 
light the bright magnitude end where the total number of stars dominates 
the source population. 



performs all galaxy fitting analysis steps from object detection to 
catalogue creation automatically. This includes (i) source detection 
and extraction with SExtractor ; (ii) preparing a ll detected objects 
for Sersic fitting with GALFIT jPeng e t al. 2002): i.e., constructing 
bad pixel masks, measuring local background levels, and setting 
up starting scripts with initial parameter estimates; (iii) running the 
Sersic model fits; and (iv) compiling all information into a final 
catalogue. 

Based on a single startup script, GALAPAGOS first runs SEx- 
tractor in the dual high dynamic range mode described in <j2.3l 
As no SExtractor setup is ever 100% optimal, we manually in- 
spected all 80 tiles for unwanted detections or over-deblended ob- 
jects. GALAPAGOS allows for the removal of such extraction fail- 
ures automatically given an input coordinate list. Additionally, we 
also composed a list of detections that are bright enough to influ- 
ence the fitting of neighbouring astronomical sources (e.g. diffrac- 
tion spikes from bright stars). Unlike the aforementioned bad de- 
tections these are not removed instantly, but kept in the source cat- 
alogue throughout the fitting process and removed only from the 
final object catalogue. Again, GALAPAGOS performs this opera- 




0.0 



15 



20 



25 



30 



[mag] 



Figure 4. Star-galaxy separation. We define a line in the magnitude-size 
plane to separate stars and galaxies (solid line). Objects above this line are 
extended galaxies; objects below are other compact objects (including most 
AGN). Grey pluses indicate all detections; black crosses only those with 
a COMBO-17 cross-match and J? a p < 24 and a redshift z > 0. Note, 
a significant number of mostly late-type stars are misidentified as galaxies 
by COMBO-17 photometry alone. The dashed line shows a line of constant 
surface brightness, which is almost parallel to our selection line at the bright 
end. 



tion automatically given a second list of coordinates. Further details 
on the process of manual fine-tuning of detection catalogues can be 
found in Barden et al. (in prep.). 

After the second run GALAPAGOS uses the cleaned output 
source list (described in $23\ to cut postage stamps for every ob- 
ject. Postage stamps are required for efficient Sersic profile fitting 
with GALFIT. The sizes of the postage stamps are based on a multi- 
ple m of the product of the SExtractor parameters KRON_RADIUS 
and AJMAGE. We define a "Kron-ellipse" with semi-major axis 
rj£ as 



r K 



m x KRON_RADIUS x AJMAGE. 



(2) 



The sky level is calculated for each source individually by 
evaluating a flux growth curve. GALAPAGOS uses the full sci- 
ence frame for this purpose in contrast to simply working on the 
postage stamp. Although in principle the background estimate pro- 
vided by SExtractor could have been used, tests show that using 
the more elabor ate GALAPAGOS scheme results in more robust 
parameter fits (Haussler et al. 2007). For a detailed description of 
the algorithm we refer to Barden et al. (in prep.). One might argue 
that GALFIT allows fitting the sky simultaneously with the science 
object. However, this requires the size of the postage stamp to be 
matched exactly to the size of the science object. If the postage 
stamp is too small, the proper sky value cannot be found; if it is too 
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V 6oe [ ma g] 

Figure 5. Fraction of extended STAGES objects and COMBO- 17 coun- 
terparts. The grey line shows the extended source fraction in STAGES. At 
bright magnitudes most sources are compact, while at the faint end almost 
all are extended. The black dotted line shows extended sources in STAGES 
with a COMBO-17 counterpart. At V 6 o6 ~ 26 the COMBO-17 complete- 
ness limit is reached. Almost no fainter sources are found in COMBO-17. 
The black solid line shows extended sources in STAGES with a COMBO- 
17 counterpart having Rap < 24. Out to Veoe ~ 22 almost every extended 
STAGES object has a COMBO-17 counterpart: the cross-correlation com- 
pleteness defined with respect to the STAGES catalogue is almost 100% 
(i.e. the ratio of black and grey lines); at Veoe ~ 24 it is ~90%. See £|3.2I 
for further discussion. 



big, computation takes unneccesarily long. Too many secondary 
sources would have to be included in the fit and the inferred sky 
value might be influenced by distant sources. Additionally, galax- 
ies may not be perfectly represented by a Sersic fit, and the sky may 
take on unrealistic values as a result. Although this method may be 
the easiest option for manual fitting, in the general case of fitting 
large numbers of sources automatically the most robust option is to 
calculate the sky value beforehand an d keep its value fix ed when 
running GALFIT (as demonstrated in Haussler et al.ll2007l) . 

Another crucial component for setting up GALFIT is deter- 
mining which companion objects should be included in the fit. 
In particular, in crowded regions with many closely neighbour- 
ing sources the fit quality of the primary galaxy improves dramati- 
cally when including simultaneously fitting Sersic models to these 
neighbours rather than simply masking them out. GALAPAGOS 
makes an educated guess as to which neighbours should be fitted 
or masked (see Barden et al., in prep, for further details). The de- 
cision is made by calculating whether the Kron-ellipses of primary 
and neighbouring source overlap. This calculation is performed not 
only for sources on the postage stamp, but on all objects on the 
science frames surrounding the current one, in order to take ob- 
jects at frame edges into account properly. Detections not identified 
as overlapping secondary sources are treated as well. Such non- 
overlapping companions are masked based on their Kron-ellipse 
and thus excluded from fitting. 

An additional requirement for fitting with GALFIT is an input 
PSF We constructed a general high S/N PSF for STAGES by com- 
bining all stars (i.e. classified by COMBO-17 photometry and hav- 
ing ACS SExtractor stellarity index > 0.85) in the brightness inter- 



val 19.5 < Vso6 ^ 23.5 and lying away from the chip edges. This 
selects non-saturated stars that can still contribute signal in their 
centres. All stars were visually inspected against binarity, compan- 
ions, or defects, which resulted in either a manually created mask, 
or the star being excluded if masking would not have been suffi- 
cient to isolate the star. With this selection 1 024 stars remained and 
were combined after subpixel cocentering and local background re- 
moval. 

In order to sample the field-variations of the PSF well and not 
be dominated by the few brightest stars, we weighted all stars iden- 
tically in the centre (where all stars carry information), but applied 
a suppression of the noise in the outer parts by a Gaussian down- 
weighting. The contribution from fainter stars in this process was 
suppressed at smaller radii relative to brighter ones. In this way we 
created a high S/N true mean PSF image of 255x255 pixel cen- 
tred exactly on the PSF and used this for all galaxy-related (but not 
AGN related) analyses. 

In its curre nt version, G ALAPAGOS sets up GALFIT to fit 
a Sersic model dSersidll968 ) for each object. A Sersic profile is 
a generalised de Vaucouleurs model with variable exponent n, the 
Sersic index: 

E (r) = E e x exp (-« [(r/r e ) 1/n - l] ) , (3) 

with the effective radius r e , the effective surface density E e , the 
surface density as a function of radius E (r) and a normalisation 
constant k = k (n). An exponential profile has n = 1 while a de 
Vaucouleurs profile has n = 4, The parameters that go into the 
model are the position [x,y], total magnitude m, the effective ra- 
dius r e , the Sersic index n, the axis ratio q (q = b/a; the ratio of 
semi-minor over semi-major half-axis ratio) and the position an- 
gle 9. Starting guesses for all parameters aside from n and r e are 
taken directly from the SExtractor output. GALAPAGOS converts 
the FLUX_RADIUS from SExtractor to estimate the effective ra- 
dius as r e = IO^'^FLUXJIADiUS 1 - 87 . This formula was found 
empiric ally to work best for simulated Sersic profiles in the GEMS 
project l lHaussler et al. 2007). The Sersic index is started at a value 
n = 2.5. 

For computational efficiency we apply constraints to the pa- 
rameter range during the fitting process. Of course, this procedure 
is not advisable when fitting objects manually, yet it is mandatory 
for an automated process like GALAPAGOS. Our constraints are 
listed in Table [3] Non-zero lower boundaries for r e and n were 
imposed for computational reasons. The maximum for r e allows 
fitting the largest galaxy in the field (750 pix correspond to ~ 60 
kpc at the cluster distance). The upper limit for the Sersic index 
is far from the de Vaucouleurs case and includes even the steepest 
profiles. The magnitude constraint flags catastrophic disagreements 
between the two photometry codes, where one of the two does not 
return a sensible result. Such problem objects may include LSB 
galaxies, where SExtractor fails to see large fractions of the total 
flux; or intrinsically faint objects with a peculiar neighbour or back- 
ground structure, where GALFIT tries to remove the excess flux. 
Objects whose values stall at the constraint limits are most likely 
not well represented by a single Sersic profile (e.g. stars or extreme 
two-component galaxies with a LSB disk). 

Finally, GALAPAGOS combines the SExtractor and GALFIT 
results into one FITS-table. At this stage flagged objects (like stel- 
lar diffraction spikes, etc.) are removed from the table. A very de- 
tailed description of GALAPAGOS including setup and computa- 
tional efficiency will be presented together with the publication of 
the code in Barden et al. (in prep.). We note that the GALFIT re- 
ported errors are purely statistical (ie. based on the assumption that 
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Table 3. GALFIT fitting constraints. 



Parameter 


Lower limit 


Upper limit 


r e 


0.3 


750 


n 


0.2 


8 


|msEx - "iGALFIT 




5 




18 20 22 24 26 
V 6oe [mag] 



Figure 6. GALFIT quality. Top panel: The two grey histograms show the to- 
tal number of fitted galaxies (light grey) and galaxies with 'bad' fits where 
the fitting procedure failed (dark grey). The heavy solid and dotted his- 
tograms show the same but for the science sample with i?ap < 24 (i.e. 
objects with a COMBO-17 counterpart only) within the overlap region 
of STAGES and COMBO-17. Bottom panel: Fraction of 'bad' fits (plus 
lcr error bars) for all fitted galaxies (dotted histogram) and those with a 
COMBO-17 match and R ap < 24 (solid histogram). Overall, ~ 23% of 
all fits ran into a constraint (dashed dotted line). For the science objects 
(STAGES/COMBO-17 cross-matched galaxies with i? a p < 24) the frac- 
tion is considerably lower (~ 8%; dashed line). The vertical line roughly in- 
dicates the surface brightness completeness limit. The 'bump' at Vgog ~ 21 
possibly results from merging two SExtractor setups (the 'hot' and 'cold' 
configurations described in £|2.31 . 



Poisson noise dominates the uncertainties of the fit parameters), 
and as such certainly under-represent the true uncertainties. A more 
meaningful measure o f uncertainties comes from fitting simulated 
galaxies, as shown in Haussler et al. I J2007I) and explored here in 
detail in ^231 

With our setup we were able to achieve an overall total of 
~ 92% high quality fits for our science targets, i.e. galaxies with 
a cross-match in the COMBO-17 catalogue and Rap < 24. We 
define 'bad' fits as those where GALFIT stalled at one of the con- 
straints in Table[3] In Fig.[6]we show the fraction of those bad fits as 
a function of SExtractor magnitude. At the bright end (Veoe < 22), 
the fraction of failures is less than 6% and rises steadily from there. 
Only when reaching the (surface brightness) completeness limit 
(roughly at Vfjo6 ~ 24 — 25) does the fraction of failed fits reach 
(and exceed) 20%. 



2.5 Completeness and Fit Quality 

To both derive completeness maps and examine fitting quality using 
GALAPAG OS, we followed a sim ilar approach as in GEMS and as 
described in lHaussler et alj d2007t) . b ut with a different, mo re real- 
istic set of simulated data. Whereas in lHaussler et al. a small 
set of only 1 600 simulated galaxies was used to find the ideal setup 
of the fitting pipeline, we have now decided on a fitting setup using 
GALAPAGOS from the start and have carried out much more in- 
tensive tests. We created entire sets of STAGES-like imaging data 
by simulating galaxies in all 80 HST/ACS tiles. Galaxies were 
simulated as single-component Sersic profiles; multi-component 
galaxies or complicated structures such as spiral arms or bars were 
not included. 

The sample of galaxies to be simulated was derived by using 
the fits of real data as described in $2.4\ From this superset, we se- 
lected a 'galaxy sample' to be simulated by excluding both stars 
and those galaxies for which the fit failed. Magnitudes and galaxy 
sizes for the simulated galaxies were chosen according to the prob- 
ability distribution of this sample. The other simulation parameters, 
(e.g. Sersic index n and axis ratio g) were then derived by choosing 
fitting values of real galaxies at approximately the same magnitude 
and size. In this way, the simulated data have parameters as close 
as possible to the real galaxy sample. 

To cover a larger number of parameter combinations, we 
slightly smoothed these values (mag by ±1 mag, log(r e ) by ±0.25 
pix, n by ±0.5 and q by ±0.2). Care was taken to make sure that q 
and n covered sensible values (0.05 < q < 1, 0.2 < n < 8). We 
also simulated galaxies two magnitudes fainter than those found 
in the real data to be able to derive completeness maps from the 
same pipeline. Twenty sets of STAGES-like data (80 tiles each) 
were simulated using this setup. In a further 50 sets, we intro- 
duced a uniform distribution of the Sersic index over the full range 
0.2 < n < 8 over all magnitudes and sizes for 5% of galaxies. 
This imposed pedestal was required in order to fill in gaps in the 
parameter space with bad number statistics or no galaxies at all, 
and was especially important for galaxies with high n-value seen 
face-on. Both position and position angle, 8, were randomly cho- 
sen for each galaxy: thus no clustering was simulated, in contrast 
to the real data. Simulating around 107 000 objects per dataset, we 
were able to derive an object density comparable to the real data 
with a mean of 60 612 galaxies found per dataset. This compares to 
75 805 galaxies in the original GALAPAGOS output from the real 
data, with ~35 000 objects in the 'galaxy' catalogue from which 
we draw the input parameters for the simulations. 

After choosing the parameters this way, we used the same sim- 
ulation script that was described in detail in lHaussler et al.1 J2007t) 
to simulate the galaxies. The images were placed in an empty image 
which was made up by empty patches of sky from the STAGES data 
to resemble the noise properties of the real data. C onvolution was 
perfor med using a STAGES PSF. In a change to the Haussl er" et al.l 
( 2007) setup, we also simulated galaxies on neighbouring tiles (or 
closely outside the data area) to realistically model effects from 
neighbouring galaxies, as well as to examine effects of combining 
the individual SExtractor catalogues within GALAPAGOS. 

By simulating fainter galaxies than are found in the real data, 
we were not only able to test the fitting quality but also the survey 
completeness. Fig. [7] shows the completeness as derived from this 
data as a function of magnitude. The left plot shows the number of 
galaxies simulated (light grey), the number of galaxies recovered 
(dark grey) and the number of galaxies with successful fit (black; 
meaning that the fit did not run into any fitting constraints). All 
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Figure 7. Completeness as a function of magnitude. Left: the number of 
simulated galaxies (light grey), recovered by SExtractor (dark grey), and 
subsequently fit successfully by GALAPAGOS (black) as a function of in- 
put magnitude. Right: Completeness functions for SExtractor (grey) and 
GALAPAGOS (black) output. One can see that GALAPAGOS returns a 
useful result in most cases. Only for relatively faint galaxies does the fit 
run into fitting constraints for a fraction of the objects. At V606 ~ 26, the 
STAGES profile fitting is therefore 80% complete. 



three histograms are normalized by the value of the bin contain- 
ing the maximum number of simulated galaxies. In total, of the 
7 497 614 galaxies simulated, 43.4% were not found in the data us- 
ing the GALAPAGOS and SExtractor setups used to analyse the 
real STAGES data. Failed objects in general were too faint to be 
detected. A further 52.5% were successfully recovered, identified 
and fitted, and 4.0% were recovered but excluded from all plots as 
the fit ran into fitting constraints. For 305 galaxies (0.004%), the fit 
crashed and did not return a result at all. 

We additionally find 51043 galaxies (0.7% of simulated 
galaxies) that could not be identified by our search algorithm, 
which looked for the closest match within 1.0". Examination of 
these galaxies shows that they are either (a) very low surface bright- 
ness galaxies for which the SExtractor positioning was not very 
secure, or (b) two neighbouring LSB galaxies that SExtractor de- 
tected as one object, also resulting in an insecure position. 

Using the whole available simulated dataset, we can derive a 
much more detailed completeness for STAGES. Magnitude alone 
is not a good estimator for completeness, as the internal light distri- 
bution has great influence on this value. More concentrated galaxy 
profiles, such as elliptical, high-n profiles, are more likely to be de- 
tected by SExtractor than disk-like low-n profiles. In addition, the 
inclination angle plays an important role. As shown in Fig. [8] we 
can divide the galaxies in different bins of n and q and for each bin 
can estimate a 2-D completeness map showing the completeness as 
a function of both magnitude and galaxy size. By looking at each 
bin one can clearly see that the completeness is indeed a function of 
magnitude as well as size. The completeness catalogue from these 
extensive simulations will be made publicly available as part of the 
STAGES data release. With the large sample and complete cover- 
age of the parameter space populated by real galaxies, one could 
make up customized completeness maps tailored to the particular 
sample in question. 

The same is true for the fitting quality. As can be seen from 
Fig-IU the fitting behaviour is a function of both surface brightness 
and Sersic index. We only show the quality as a function of Sersic 
index, but again one can determine fitting quality as a function of 
any combination of the fitted parameters. One can see that high-n 
galaxies are harder to fit than low-n galaxies, e.g. the magnitude 
deviation A is 0.00 (a = 0.07) at around the sky level for galax- 
ies with < n < 1.5, while A = 0.03 (<r = 0.12) at the high- 



Table 5. COMBO-17 observing runs with A901/2 imaging. 



COMBO-17 ran code 


Dates 




A 


11.02. 


-22.02.1999 


E 


28.01. 


-11.02.2000 


G 


19.01. 


-20.01.2001 



est n-bin. The effect is even larger at fainter galaxies: A = 0.00 
(a = 0.18) at 25 mag arcsec~ 2 and A = 0.08 (a = 0.28) for 
low- and high-n galaxies, respectively. A similar trend can be seen 
for galaxy sizes: A = 0.8% {a = 7.3%) and A = -3.7% 
(a = 19.0%) at the sky level, and A = -0.4% (a = 18.3%) 
and A = —10.7% (a = 36.1%) at 25 mag arcsec -2 . If one exam- 
ines relative deviations of the Sersic index, there is essentially no 
trend seen between different bins of n. In an absolute sense, then, 
the Sersic index is still less well recovered in the high-n bin. 

In general, the systematic deviations are very small except at 
the faintest galaxies detectable, and both deviation A and a of the 
distrib utions are well under stood within STAGES. As was pointed 
out in lHaussler et alj (l2007). the uncertainties returned by GALFIT 
(and therefore GALAPAGOS) underestimate the true uncertainty 
by a large amount. Using a statistical approach therefore returns 
more reliable errorbars for the individual parameters. The simula- 
tions and catalogue presented here allow a flexible means of esti- 
mating errors on profile fitting for any possible subsample of galax- 
ies. 



3 COMBO-17 DATA 

3.1 COMBO-17 observations and catalogue 

In this section we briefly describe the COMBO-17 data on the 
A901/2 field, including observations, catalogue entries and object 
sa mples. The corres ponding data on the CDFS field were published 
in lWolf et al.l J20041 hereafter W04), where further technical details 
can be found. 

The filter set (Table |4) contains five broad-band filters (UB- 
VRI) and 12 medium-band filters covering wavelengths from 350 
to 930 nm. All observations were obtained with the Wide Field 
Imager (WFI) at the MPG/ESO 2.2m-telescope on La Silla, Chile. 
A field of view of 34' x 33' (see Fig. Hi is covered by a CCD 
mosaic consisting of eight 2k x 4k CCDs with a scale of O'.^SS 
per pixel. The observations on the A901/2 field were spread out 
over three observing runs between January 1999 and February 
2001. They encompass a total exposure time of ~185 ks of which 
~20 ks were taken in the .R-band during the best seeing condi- 
tions. A dither pattern with at least ten telescope pointings spread 
by Aa, AS < ±72" allowed us to cover the sky area in the gaps 
of the CCD mosaic. 

Flux calibration was done with our own tertiary standard stars 
based on spectrophotometry observations, a suitable method to 
achieve a homogeneous photometric calibration for all 17 WFI fil- 
ter bands. Two G stars with B ~ 15 (with COMBO-17 identifi- 
cation numbers 45811 and 46757) were observed at La Silla with 
DFOSC at the Danish 1.54 m telescope. A wide (5") slit was used 
for the COMBO-17 standards as well as for an external calibrator 
star. 

The object searc h for the COMBO-17 s ample was done with 
SExtractor software (Bertin & Arnouts 1996) in default setup, ex- 
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Figure 8. Completeness maps as a function of Sersic index n and axis ratio q (as labelled above and to the right of the plots). To guide the eye, we overplot 
a vertical line at mag 26 and a surface brightness line (diagonal, dashed) at 28 mag arcsec~ 2 . As one can clearly see, the completeness (shown in greyscale, 
black is complete, white is incomplete or no data) is a strong function of all magnitude, size (and therefore surface brightness), q and n. The outline contour 
shows the region in this plot where galaxies have been simulated to demonstrate where these plots are reliable. 



cept for choosing a minimum of 12 significant pixels required for 
the detection of an object. We first search rather deep and then 
clean the list of extracted objects of those having a S/N ratio below 
4, which corresponds to > 07 1 2422 error in the total magnitude 
MAG_BEST. As a result we obtained a catalogue of 63 776 objects 
with positions, morphology, total i?-band magnitude and its error. 
The astrometric accuracy is better than 0'.' 15. Using our own aper- 
ture photometry we reach a 5a point source limit of R « 25.7. 

We obtained spectral energy distributions of all objects from 
photometry in all 17 passbands by projecting the known object co- 
ordinates into the frames of reference of each single exposure and 
measuring the object fluxes at the given locations. In order to op- 
timize the signal-to-noise ratio, we measure the spectral shape in 
the high surface brightness regions of the objects and ignore po- 
tential low surface brightness features at large distance from the 
centre. However, this implies that for large galaxies at low redshifts 
2 < 0.2 we measure the SED of the central region and ignore 
colour gradients. 

Also, we suppressed the propagation of variations in the see- 



ing into the photometry by making sure that we always probe the 
same physical footprint outside the atmosphere of any object in 
all bands irrespective of the PSF. Here, the footprint f(x, y) is 
the convolution of the PSF p(x, y) with the aperture weighting 
function a(x,y). If all three are Gaussians, an identical physi- 
cal footprint can be probed even when the PSF changes, simply 
by adjusting the weighting function a(x, y). We chose to measure 
fluxes on a footprint of l'.'5 FWHM outside the atmosphere (~ 4.2 
kpc at z ~ 0.165). In detail , we use the package MPIAPHOT 
( Meisenhei mer & R6ser|[l993r) to measure the PSF on each indi- 
vidual frame, choose the weighting function needed to conserve the 
footprint and obtain the flux on the footprint. Fluxes from individ- 
ual frames are averaged for each object and the flux error is derived 
from the scatter. Thus, it takes not only photon noise into account, 
but also suboptimal flatfielding and uncorrected CCD artifacts. 

All fluxes are finally calibrated by the tertiary standards 
in our field. The aperture fluxes correspond to total fluxes for 
point sources, but underestimate them for extended sources. The 
difference between the total (SExtractor-based) and the aperture 
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Figure 9. Fit Quality. The deviations of the most important galaxy parameters as a function of surface brightness. Top row: Magnitude deviation (fit - simulated), 
middle row: Size ratio (fit/simulated), bottom row: Sersic index n (fit/simulated). Contours show the data normalized by the number of galaxies in each surface 
brightness bin; the black solid line shows the mean of the distribution; and the black dashed lines show the sigma of the distribution (3<r in case of magnitudes, 
lcr in size and Sersic index). All plots are shown for different Sersic indices as labelled above the plots. The vertical grey line represents the mean brightness 
of the sky background in STAGES. The magnitude and sizes are less well recovered in high-n galaxies, but the relative recovery of n is similar in all cases. 

Table 4. COMBO-17 imaging data on the A901/2 field: For all filters we list total exposure time, average PSF among individual frames, the 10<r (Vega) 
magnitude limits for point sources and the observing runs (see Tab. |5J in which the exposure was collected. For flux and magnitude conversions we list AB 
magnitudes and photon fluxes of Vega in all filters. The i?-band observations were taken in the best seeing conditions. 



Accn/fwhm 


^exp 


seeing 


»"lim,10<r 


run code 


mag of Vega 


Fphot of Vega 


(run) 




(sec) 




(Vega mags) 




(AB mags) 


(10 s photons/m 2 /nm/s) 


365/36 


U 


22100 


l'.'lO 


23.7 


G 


+0.77 


0.737 


458/97 


B 


20500 


l'.'20 


25.4 


A, G 


-0.13 


1.371 


538/89 


V 


6000 


1'.'20 


24.3 


E 


-0.02 


1.055 


648/160 


R 


20300 


tf'75 


25.0 


E 


+0.19 


0.725 


857/147 


I 


7500 


I'.'OO 


22.7 


E 


+0.49 


0.412 


418/27 




7300 


l'/20 


24.0 


E 


-0.19 


1.571 


462/13 




10000 


l'.'lO 


23.7 


E 


-0.18 


1.412 


486/31 




5500 


1'.'15 


24.0 


E 


-0.06 


1.207 


519/16 




6000 


l'.'OS 


23.6 


E 


-0.06 


1.125 


572/25 




5000 


<y.'85 


23.5 


E 


+0.04 


0.932 


605/21 




6000 


tf'95 


23.4 


E 


+0.10 


0.832 


645/30 




4950 


1C30 


22.7 


E 


+0.22 


0.703 


696/21 




6600 


I'.'OO 


22.7 


E 


+0.27 


0.621 


753/18 




7000 


l'.'OS 


22.2 


E 


+0.36 


0.525 


816/21 




19200 


0'.'85 


22.8 


A 


+0.45 


0.442 


857/15 




16600 


ms 


21.7 


E 


+0.56 


0.386 


914/26 




15700 


tf'95 


21.9 


E 


+0.50 


0.380 
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(MPIAPHOT-based) magnitude is listed as an aperture correction 
and used to calculate e.g. luminosities. For further details on the 
observations and the data processing, see W04. 

The A901/2 field is affected by substantial foreground dust 
reddening at the level of E(B — V) w 0.06, in contrast to the 
CDFS. Hence, any SED fitting and derivation of luminosities re- 
quires dereddened SEDs. Therefore, in the catalogue we list three 
sets of photometry: 

(i) i?-band total and aperture magnitudes as observed for the 
definition of samples and completeness; 

(ii) aperture fluxes _F p hot in 17 bands, dereddened using Ay = 
0.18 and (A v , A B ,A R , Ai) = A v X (1.63, 1.24, 0.82, 0.6) with 
similar numbers for medium-band filters (rereddening with these 
numbers would restore original measurements); and 

(iii) aperture m agnitudes (Vega) in all 17 bands, dereddened, on 
the Asinh system dLupton et al.ll9 99) that can be used for logarith- 
mic flux plots with no trouble arising from formally negative flux 
measurements. 

Fluxes are given as photon fluxes F p hot in units of 
photons/m 2 /s/nm, which are related to other flux definitions by 

vF v = hcF phot = XF\ . (4) 

Photon fluxes are practical units at the depth of current surveys. 
A magnitude of V = 20 corresponds to 1 photon/m 2 /s/nm in all 
systems (AB, Vega, ST), provided V is centred on 548 nm. Flux 
values of an object are missing in those bands where every exposure 
was saturated. 

The final catalogue contains quality flags for all objects in an 
integer column ('phot.flag'), holding the original SExtractor flags 
in bit to 7, corresponding to values from to 128, as well as some 
COMBO- 17 quality control flags in bits 9 to 11 (values from 512 
to 2048). We generally recommend that users ignore objects with 
flag values phot_flag 8 for any statistical analysis of the object 
population. If an object of particular interest shows bad flags, it 
may still have accurate COMBO- 17 photometry and could be used 
for some purposes. Often only the total magnitude was affected by 
bright neighbours, while the aperture SED is valid. 

We then employ the usual COMBO- 17 classification and red- 
shift estimation by template fitting to libraries of stars, galaxies, 
QSOs and white dwarfs. There, the error rate increases very signif- 
icantly at i? ap > 24. We refer again to W04 for details of the li- 
braries and known deficiencies of the process, but repeat here (and 
correct a misprint in W04) the definition of the classifications (see 
Table|6j. 

We also show in Table [6] a comparison of the sample sizes 
in different classes between the A901/2 and the CDFS field of 
COMBO-17. The main difference is that the A901/2 field contains 
more than twice the number of stars given its position at relatively 
low galactic latitude (+33.6 deg). Another difference is that it con- 
tains 30% more galaxies than the CDFS, which is both a conse- 
quence of the cluster A901/2 and the underdensity in the CDFS 
seen at z ~ [0.2, 0.4]. Fig.[l0]shows a colour-magnitude diagram 
of the star and white dwarf sample as well as redshift-magnitude 
diagrams for galaxies and QSOs. 

Redshifts are given as Maximum-Likelihood values (the peak 
of the PDF), or as Minimum-Error- Variance values (the expectation 
value of the PDF). MEV redshifts have smaller true errors, but are 
only given when the width of the PDF is lower than a z /(l + z) < 
0.125. If PDFs are bimodal with modes of sufficiently small width, 
then both values are given with the preferred (larger-integral) mode 



Table 6. Definition of entries for the 'mc.class' column and comparison of 
object numbers between the COMBO-17 data sets of the A901/2 and CDFS 
field. The samples refer to a magnitude range of i? ap = [16, 24] and only 
objects with phot_flag< 8. The A901/2 field is richer in stars because of its 
galactic coordinates. It is also richer in galaxies due to the cluster, while the 
CDFS is underdense at z = [0.2,0.4]. We note that these definitions are 
based on the COMBO-17 data SED and morphology; star-galaxy separation 
employing morphological information from the HST imaging (Equation fT} 
is considered separately. 



Class entry 


Meaning 


ALA901/2 


ALCDFS 


Star 


stars 

(only point sources) 


2096 


992 


WDwarf 


white dwarf 

(only point sources) 


14 


9 


Galaxy 


galaxies 

(shape irrelevant) 


14555 


11054 


Galaxy (Star?) 


binary or low-z galaxy 
(star SED but extended; 
ambiguous colour space) 


44 


46 


Galaxy (Unci!) 


SED fit undecided 
(most often galaxy) 


316 


243 


QSO 


QSOs 

(only point sources) 


73 


66 


QSO (Gal?) 


Seyfert-1 AGN or 
interloping galaxy 
(AGN SED but extended; 
ambiguous colour space) 


36 


31 


Strange Object 


unusual strange spectrum 

(X? cd > 30) 


1 


3 



providing the primary redshift. Our team uses only MEV redshifts 
(with column name 'mc_z') for their analyses. 

The galaxy sample with MEV redshifts is > 90% complete 
at all redshifts for i? ap < 23. Near z ~ 1, the MEV redshifts are 
this complete even at i? ap = 24. Below this cut, increasing photon 
noise drives an expansion of the width of the PDF. The error limit 
for MEV redshifts then makes the completeness of galaxy sam- 
ples with MEV redshifts drop. The 50% completeness is reached 
at ii ~ 24 to 25 depending on redshift. These results have been 
determined from simulations and are detailed in W04. Complete- 
ness maps are included in the data release and take the form of a 
3-D map of completeness depending on aperture magnitude, red- 
shift and restframe U — V colour. 

To date, the photo-z quality on the A901/2 field has only been 
investigated with a comparison to spectroscopic redshifts at the 
bright end. W04 reported results from a sample of 404 bright galax- 
ies with R < 20 and z = [0, 0.3], 351 of which were on the A901/2 
field, and 249 of which were members of the A901/2 cluster com- 
plex (t) |4.5t . The other 53 objects were observ ed by the 2dFGRS on 
the CDFS and SI 1 fields JColless et al.l200lh . There we found that 
77% of the sample had photo-z deviations from the true redshift 
\S Z /(1 + z)\ < 0.01. Three objects (less than 1%) deviate by more 
than 0.04 from the true redshift. 

Currently, we do not have faint spectroscopic samples on the 
A901/2 field, however a spectroscopic dataset from VVDS exists 
on the COMBO-17 CDFS field. From a sample of 420 high-quality 
redshifts that are reasonably complete to i? ap < 23, we find a la 
scatter in S z /(1 + z) of 0.018, but also a mean bias of —0.011. Fur- 
thermore, the faint CDFS data show ~ 5% o utliers with deviations 
of more than 0.06 (Hildeb randt et alj|2008h . From a collection of 
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Figure 10. Left panel: Stars (dots) and white dwarfs (crosses): B — V colour vs. i?tot- The two reddest stars at R ss 23 and B — V > 2 are M5-6 stars. 
Centre panel: Red-sequence (black) and blue-cloud galaxies (green): MEV redshift vs. i?tot- Right panel: QSOs: MEV redshift vs. iJtot- 



Table 7. The restframe passbands and their characteristics. 



name 


Acen/fwhm 


mag of Vega 


Fphot of Vega 




(nm) 


(AB mags) 


(10 8 phot/m 2 /nm/s) 


(synthetic) 


145/10 


+2.33 


0.447 


(synthetic) 


280/40 


+1.43 


0.529 


Johnson U 


365/52 


+0.65 


0.820 


Johnson B 


445/101 


-0.13 


1.407 


Johnson V 


550/83 


+0.00 


1.012 


SDSSu 


358/56 


+0.84 


0.704 


SDSS g 


473/127 


-0.11 


1.305 


SDSSr 


620/115 


+0.14 


0.787 



spectroscopic samples we modelled the overall ler redshift errors 
at R < 24 and z < 1 in W04 as 

a z /(l + z) M 0.005 x v 7 ! + lO - 6 ^- 20 - 5 ) . (5) 

Later we use a variant of this approximation to estimate the com- 
pleteness of photo-z based selection rules for cluster members. 

The template fitting for galaxies produces three parameters, 
i.e. redshift as well as formal stellar age and dust reddening val- 
ues. The age is encoded in a template number running from 
(youngest) to 59 (oldest), where we use the same PEGASE (see 
iFioc & Rocca-Vo lmerange 1997, for discussion of an earlier ver- 
sion of the model) template grid as described in W04. The look 
back times to the onset of the r = 1 Gyr exponential burst range 
from 50 Myr to 15 Gyr. 

Restframe properties are derived for all galaxies and QSOs as 
described in W04. Table [7] lists the restframe passbands we cal- 
culate and gives conversion factors from Vega magnitudes to AB 
magnitudes and to photon fluxes. The SED shape is defined by the 
aperture photometry and the overall normalization is given by the 
total SExtractor photometry from the deep i?-band. However, if a 
galaxy has both a steep colour gradient and a large aperture correc- 
tion, then the restframe colours will be biased by the nuclear SED. 

The column 'ApD_Rmag' contains the magnitude difference 
between the total object photometry and the point-source cali- 
brated, seeing-adaptive aperture photometry: 



ApD_Rmag — Rmag — Ap_Rmag . (6) 

On average, this value is by calibration zero for point sources, and 
becomes more negative for more extended sources. 

3.2 Cross-correlation of STAGES and COMBO-17 
catalogues 

Having created separate catalogues from the STAGES ( i]2.3| i ]2.4b 
and COMBO-17 ( i]3.U datasets, we next wish to create a com- 
bined, master catalogue. In GEMS, this was accomplished by ap- 
plying a nearest neighbour matching algorithm with a maximum 
matching radius of 0775. The choice of maximum radius is gov- 
erned by the resolution of the two datasets (HST: 0" 1 ; COMBO-17: 
0775). 

For STAGES we have however chosen to improve over this 
approach. For most galaxies, their measured centres do not change 
if the input image is smoothed. For example, if the HST image 
of a normal spiral or elliptical galaxy is convolved with a Gaus- 
sian function to match the ground-based seeing, the centre esti- 
mated from the high-resolution (in this case STAGES) and the 
low-resolution (here COMBO-17) images should coincide. For dis- 
torted galaxies or mergers, this may no longer be the case. Instead, 
the brightest peak in the STAGES image, detected as the object 
centre by SExtractor, may be relatively far from the centre in the 
COMBO-17 image. 

In order to maximise the number of good matches between 
STAGES and COMBO-17, in particular at low redshift, i.e. A901/2 
cluster distance, we have devised the following scheme. For 
STAGES the average source density corresponds to roughly two 
objects per 5"-radius circle. We cross-correlate the STAGES and 
COMBO-17 catalogues using a nearest neighbour matching algo- 
rithm as described above with a maximum matching radius of 5". 
The resulting matches we plot in Fig. QT] (left panel). In particular 
at faint magnitudes many matches are found that appear unrelated. 
In contrast, at brighter magnitudes several sources are correlated at 
radii much larger than the COMBO-17 seeing (0.75"), which still 
identify the same object. In Fig.QT]we also show a line that subdi- 
vides the plot into two regions: 

d m = -0.3 x (Veoe - 29) , (7) 

with the matching radius d m in arcsec and the STAGES SExtractor 
magnitude Veo6- Below the line, objects are considered to be cor- 
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Figure 11. Cross-correlation of HST and COMBO-17 data. Left: The distance to the nearest neighbour within a search radius of 5"is plotted as a function 
of HST magnitude. At the faint end galaxies are matched to uncorrelated neighbours. Resolving irregular structures in the HST images results in detected 
galaxy centres being located farther from the COMBO-17 galaxy centre than a seeing distance. Matching bright objects at large separations while removing 
random correlations at faint fluxes requires a cut as indicated by the diagonal line. Objects within the box (V606 < 25 and 1" < match distance < 2.5") were 
inspected by eye. Right: Ratio of matching distance and Kron size as a function of HST magnitude. Values larger than ~ 1 imply a matching radius larger 
than the object size in the HST image. Sources with Rap < 24 are shown as black symbols; objects with a match below the cut (diagonal line in left panel) 
are plotted in dark grey; the remaining sources with a match within 5" are shown as light grey symbols. 



8000 




Angular separation [arcsec] 



Figure 12. Histogram of matching radii for all objects (outer histogram) 
and i? ap < 24 objects (inner histogram). The typical angular separation 
between a COMBO-17 object with i? a p < 24 and its HST counterpart is 
~ 0.12" ±0.08". 



is uncorrelated. From the STAGES SExtractor data we estimate the 
'extent' of an object by its Kron size K = rj^ x a, from the Kron 
radius and semi-major axis radius a. We limit the Kron size 
to K > 0.75". A ratio of d m /K > 1 indicates that the matched 
COMBO-17 source lies outside the region covered by the object 
in the STAGES image. In Fig. [TT] (right panel) we overplot in grey 
all sources that were assigned a partner from the nearest neighbour 
matching. This provides further evidence for the improved quality 
of our new cross-correlation method. 

In summary, the combined catalogue contains 88 879 sources. 
Of these, ~ 6 577 objects with a COMBO-17 ID are not within 
the region covered by the STAGES HST mosaic (~ 1 664 of 
these have 7? a p < 24). Moreover, ~ 1 271 STAGES detections 
are outside the COMBO-17 observation footprint[f| Inside the re- 
gion covered by both surveys, there are ~ 81 031 sources. For 
50701 objects the method described above provides a match be- 
tween COMBO-17 and STAGES (15760 of these have i?ap < 24). 
~ 23 833 sources detected in STAGES do not have counterparts 
in COMBO-17; ~ 6 497 sources from the COMBO-17 catalogue 
are not matched to STAGES detections. Out of these, only ~ 79 
objects have i?ap < 24. We therefore emphasize that for our sci- 
ence sample of COMBO-17 objects, defined as having Rap < 24, 
99.9% have a STAGES counterpart. The majority of failures result 
from confusion by neighbouring objects or simply non-detections. 



related, while above they are not correlated. This division is em- 
pirically motivated by the requirement to match objects at the faint 
end out to the COMBO-17 resolution limit (0.5"-1.0") while also 
correlating sources at larger radii at the bright end. The slope of 
the curve was determined by visual inspection of the matches in- 
side the indicated box. Typically, the distance between centroids is 
~tf'l(Fig.[2). 

Another way of investigating this issue is by calculating 
whether the nearest matching neighbour falls within the area cov- 
ered by the object in the STAGES image. If the projected COMBO- 
17 position is beyond the optical extent of the source in STAGES, it 



3.3 Selection of an A901/2 cluster sample 

We wish to define a 'cluster' galaxy sample of galaxies belonging 
to the A901/2 complex for various follow-up studies of our team 
that are in progress. These studies may have different requirements 

3 The observation footprint for both STAGES and COMBO-17 is rather 
difficult to determine. Therefore, we provide only approximate numbers 
good to ~ 50 objects. A more elaborate scheme than the one used to pro- 
duce these numbers is well beyond the scope of this paper. 
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for the completeness of cluster members and the contamination by 
field galaxies. We therefore quantified how these two key values 
vary with both magnitude and width of the redshift interval in order 
to inform our choice of definition. 

The photo-z distribution of cluster galaxies was assumed to 
follow a Gaussian with a width given by the photo-z scatter in 
Equation [5] The distribution of field galaxies was assumed to be 
consistent with the average galaxy counts n(z, R) outside the clus- 
ter and varies smoothly with redshift and magnitude assuming no 
structure in the field. Samples were then defined by redshift inter- 
vals Zphot = [0.17 — Az, 0.17+ Az], where the half-width Az was 
allowed to vary with the magnitude^ We calculated completeness 
and contamination at all magnitude points simply using the counts 
of our smooth models. 

We found that as long as the half-width in redshift is not 
much larger than a couple of Gaussian FWHMs, the contamina- 
tion changes only little. The ratio of selected cluster to field galax- 
ies is almost invariant as shrinking widths cut into numbers for 
both origins. Only enlarging the width significantly over that of the 
Gaussian increases contamination by field galaxies. On the con- 
trary, such large widths do not affect the completeness of the clus- 
ter sample much, while shrinking the width too far eats into the true 
cluster distribution and reduces completeness of the cluster sample. 

For our purposes, we compromised on a photo-z width such 
that the completeness is > 90% at any magnitude, just before fur- 
ther widening starts to increase the contamination above its mag- 
dependent minimum (see Fig.[T3]and Fig. [14] left panel). For this 
we chose a half-width of 

Az(R) = ^0.015 2 +0.0096525 2 (1 + io°-6(«tot-20.s)) . (g) 

This equation defines a half-width that is limited to 0.015 at 
the bright end and expands as a constant multiple of the estimated 
photo-z error at the faint end. The floor of the half-width is mo- 
tivated by including the entire cluster member sample previously 
studied by WMG05. The completeness of this selection converges 
to nearly 100% for bright galaxies, as a result of intentionally in- 
cluding the WGM05 sample entirely. 

The right panel of Fig. [14] shows that the differential contam- 
ination increases rapidly towards faint magnitudes, simply as a re- 
sult of the photo-z error-driven dilution of the cluster sample. Here, 
contamination means the fraction of galaxies that are field mem- 
bers, as measured in a bin centred on the given magnitude with 
width 0. 1 mag. Contamination at a given apparent magnitude trans- 
lates into contamination at a resulting luminosity at the cluster dis- 
tance (except that scatter in the aperture correction smears out the 
contamination relation slightly). 

Already at R ap = 23.2 the sample contains as many cluster 
as field members. This corresponds to My ~ —16.5 for the aver- 
age galaxy, but scatters around that due to aperture corrections. As 
we probe fainter this selection adds more field galaxies than clus- 
ter members. Follow-up studies can now determine an individual 
magnitude or luminosity limit given their maximum tolerance for 
field contamination. For example, WGM05 selected cluster galax- 
ies at My — < —17.775 (My < —17 for their adopted cosmology 
with Ho = 100 km s" 1 Mpc" 1 ) for an earlier study of the A901/2 
system in order to keep the contamination at the faint end below 
20%. 



4 We use .Zphot = 0.17 for the mean cluster redshift here rather than the 
spectroscopically confirmed z spC c ~ 0.165 due to the known bias dis- 
cussed in £|3.1I 
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Figure 13. MEV redshift estimate vs. total i?-band magnitude. The 'galaxy' 
sample is shown in green, while the sample of 'cluster' galaxies defined by 
Equation [8] is shown in black. The magnitude-dependent redshift interval 
guarantees almost constant high completeness, while the field contamina- 
tion increases towards faint levels (Fig. 1141 . We note that at faint magnitudes 
there is an apparent asymmetry towards lower redshift at faint magnitudes 
within the cluster sample. The photometric redshifts may be skewed by sys- 
tematic effects but the average —0.02 offset at i?tot ~ 22.5 is within the 
lcr error envelope. 
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Figure 14. Left panel: Completeness of the cluster sample defined in Fig. l 131 
and designed to provide high completeness at all magnitudes. Right panel: 
The field contamination of the cluster sample increases at faint levels due 
to photo-z dilution of the cluster. Narrowing the selected redshift interval 
would not reduce the contamination. Contamination rates are estimated to 
be (10, 20, 30, 50, 70)% at R ap = (20.3, 21.65, 22.3, 23.2, 24.0). 



The cluster sample thus obtained covers quite a range of 
photo-z values at the faint end, and restframe properties are derived 
assuming these redshifts to be correct. However, if we assume a 
priori that an object is at the redshift of the cluster, then we may 
want to know these properties assuming a fixed cluster redshift of 
2 = 0.167. Hence, the SED fits and restframe luminosities are re- 
calculated for this redshift and reported in additional columns of 
the STAGES catalogue in Table IB 1 1 (with '_cl' suffix indicating 
cluster redshift). Of course, if the a-priori assumption is to believe 
the redshifts as derived, then the original set of columns for which 
we have derived the values is relevant. 
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4 FURTHER MULTIWAVELENGTH DATA AND 
DERIVED QUANTITIES 

In this section we describe further multiwavelength data for the 
A901/2 region taken with other facilities (Fig. [2}. We also present 
several resulting derived quantities (stellar masses and star forma- 
tion rates) that appear as entries in the STAGES master catalogue. 

4.1 Spitzer 

Spitzer observed a 1° x 0?5 field around the A901/2 system in 
December 2004 and June 2005 as part of Spitzer GO-3294 (PI: 
Bell). The MIPS 24u.m data were taken in slow scan-map mode, 
with individual exposures of 10 s. We reduced the individual image 
frames using a custom dat a-analysis tool (DAT) developed by the 
GTOs iGordon et al.l2005h . The reduced images were corrected for 
geometric distortion and combined to form full mosaics; the reduc- 
tion which we currently use does not mask out asteroids and other 
transients in the mosaicing0 The final mosaic has a pixel scale of 
1'.'25 pixeP 1 and an image PSF FWHM of ~ 6". Source detec- 
tion and photometry w ere performed using techniques described in 
IPapovich et all |2004); based on the analysis in that work, we es- 
timate that our source detection is 80% complete at 97 ^JjQ for a 
total exposure of ~ 1400 spix -1 . By detecting artificially-inserted 
sources in the A901 24 image, we estimated the completeness of 
the A901 24 fim catalog. The completeness is 80%, 50% and 30% 
at 5, 4 and 3a, respectively. 

Note that there is a very bright star at 24p.m near 
the centre of the field at coordinates (a, <5)j20oo = 
(09 ,l 56 m 32 i :4, -10°01'15") (see SjAl] for details of this ob- 
ject). In our analysis of the 24u.m data we discard all detections 
less than 4' from this position in order to minimise contamination 
from spurious detections and problems with the background level 
in the wings of this bright star. It is to be noted that there are a 
number of spurious detections in the wings of the very brightest 
sources; while we endeavoured to minimise the incidence of these 
sources, they are difficult to completely eradicate without losing 
substantial numbers of real sources at the flux limit of the data. 

To interpret the observed 24(j.m emission, we must match the 
24p.m sources to galaxies for which we have redshift estimates 
from COMBO-17. We adopt a l"matching radius. In the areas 
of the A901/2 field where there is overlap between the COMBO- 
17 redshift data and the full-depth MIPS mosaic, there are a total 
of 3506(5545) 24u.m sources with fluxes in excess of 97(58)/iJy. 
Roughly 62% of the 24u.m sources with fluxes > 58/iJy are de- 
tected by COMBO-17 in at least the deep i?-band, with R < 26. 
Some 50% of the 24(im sources have bright R to t < 24 and have 
photometric redshift z < 1; these 50% of sources contain nearly 
60% of the total 24(j.m flux in objects brighter than 58/xJy. Sources 
fainter than R > 24 contain the rest of the /24 > 58/iJy 24(j.m 
sources; investigation of COMBO-17 lower confidence photomet- 
ric redshifts, their optical colours, and results from other studies 
lends weight to the argument that essentially all of these sources 
are at z > 0.8, with the bu lk l ying at z > 1 (e.g.lLe Floc'h et al.l 
|2004 IPapovich et alj|2004l ; see|Le Floc'h et al. 2005 for a further 

5 This only minimally affects our analyses because we match the IR de- 
tections to optical positions, and most of the bright asteroids are outside the 
COMBO-17 field. 

6 We note that for previous papers we used the catalogue to lower flux 
limits, down to 3<r; accordingly, we have included such lower-significance 
(and more contaminated) matches in the catalogue. 



discussion of the completeness of redshift information in the CDFS 
COMBO-17 data). 

Observations with IRAC (Infrared Array Camera; iFazio et all 

120041) at 3.6, 4.5, 5.8 and 8.0u.m were also taken as part of this 
Spitzer campaign: those data are not discussed further here, and 
will be described in full in a future publication. 

4.2 Star formation rates 

We provide estimates of star formation rate, determined using 
a combination of 24[i.m data (to probe the obscured star for- 
mation) and COMBO-17 derived rest-frame 2800A luminosities 
(to probe unobscured star formation). Ideally, we would have a 
measure of the total thermal IR flux from 8-1000u.m; instead, 
we have an estimate of IR luminosity at one wavelength, 24u.m, 
corresponding to rest-frame 22-1 2 u.m at the redshifts of interest 
2 = 0.1 — 1. Local IR-luminous galaxies show a tight correla- 
tion bet ween rest-frame 12-15[i.m luminosity and total IR luminos- 
ity (e . g-ISpinoglio et alll995t ICharv & Elbazll200ll : iRoussel et ail 
l200lt|papovich & Bellll2002l), w ith a scatter of ~ 0.15 dexQ Fol- 
lowing Pa povich & Bell {2002), we choose to construct total IR 
luminosity from th e observed-fram e 24u,m data. We use the Sbc 
template from the iDevriendt et al . (1999) SED library to trans- 
late observed-frame 24u.m flux into the 8-lOOOp.m total IR lu- 
minosity[f] The IR luminosity uncertainties are primarily system- 
atic. Firstly, there is a natural diversity of IR spectral shapes at a 
given galaxy IR luminosity, stellar mass, etc.; one can crudely esti- 
mate the sc ale of this uncer t ainty b y using the full range of tem- 
plat es from Devriendt et alj dl999l) . or by using templates from, 
e.g. jDale et alj (2001) instead. This uncertainty is <0.3dex (this 
agrees roug hly with the scatter s een between 24u.m luminosity and 
SFR seen in Calz etti et al .l2007h . Secondly, it is possible that a sig- 
nificant fraction of 0.1 < z < 1.0 galaxies have IR spectral energy 
distributions not represented in the local Universe: while it is im- 
possible to quantify this error until the advent of Herschel Space 
Telescope, current results suggest that the bulk of intermediate- 
high redshift gal axies have IR spec t ra similar to galaxies in the 
local universe jAppleton et al.ll2004l : lElbaz et alj|2005t lYan et alj 
l2005l : IZheng etal] 20071) . 

We estimate SFRs using the combined directly-observed UV 
light from young stars and the dust-reproc essed IR e mission of 
the sam ple galaxies (e.g. jGordon et al.l2 000). Following Be ll et all 
(2005), we estimate SFR ip using a calibration derived from PE- 
GASE ass uming a 100 M yr-old stellar population with constant 
SFR and a lChabried < |2003|) IMF: 

t/7 (M yr" 1 ) = 9.8 x 10" 11 x (L m + 2.2L UV ), (9) 

where Lxr is the total IR luminosity (as estimated above) and 
Luv = 1-5i/Zi/,2800 is a rough estimate of the total integrated 
1216A-3000A UV luminosity, derived using the 2800A rest-frame 
luminosity from COMBO-17 L,2soo- The factor of 1.5 in the 
2800A-to-total UV conversion accounts for the UV spectral shape 
of a 100 Myr-old population with constant SFR, and the UV flux 



7 Star-forming regions in local galaxies appear to follow a slightly non- 
linear re lation between rest- frame 24irm emission and SFR, with SFR oc 
^24^im ^Calzetti et alj2007h . although note that this calibration is between 
24irm emission and SFR (not total IR luminosity). 

8 Total 8-1000inn IR luminosit ies are ~ . 3 dex higher than the 42.5- 
122.5 |^m luminosities defined by H elou et all J 1 988t> . with an obvious dust 
temperature dependence. 
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is multiplied by a factor of 2.2 before being added to the IR lumi- 
nosity to account for light emitted longwards of 3000A and short- 
wards of 1216A by the unobscured young st ars. This SFR c alibra- 
tion is derived using identical assumptions to lKennicuttl J 1998h . and 
the calibration is consistent with his to within 30% once different 
IMFs are accounted for. Uncertainties in these SFR estimates are 
a factor of two or more in a galaxy-by-galaxy sense, and system- 
atic uncertainty in the overall SFR scale is l ikely to be less than a 
factor of two (see, e.g.. lBelllp2003l ; lBell et ai]|2005l for further dis- 
cussion of uncertainties). The adopted calibration assumes that the 
infrared luminosity traces the emission from young stars only; con- 
tributions from potential AGN can be identified an d excluded by 
cross- matc hing with the X-ray a nd optical data as in lGilmour et al.1 
d2007h and lGallazzi et alj J2008h . 

Again, for galaxies in the 'cluster' sample, we present also 
SFR estimates assuming that the galaxies are at the cluster redshift 
with the suffix '_cl' in added to the column name. 



4.3 Stellar Masses 



Borch et al. I d2006h estimated the stellar masses of galaxies in 
COMBO- 17 using the 17-passband photometry in conjunction 
with a template library derived using the PEGASE stellar popu- 
lation model. The non-evolving template stellar populations had an 
age/metallicity combination equivalent to roughl y solar metallicity 
and ~ 6 Gyr since the s t art of star formation j B orch et all 1 2006) 
adopte d a Kroupa et alj d 19931) stellar IMF; the use of a [kroupa 
d200lh or lChabrieJd2003h IMF would have yielded the same stellar 
masses to within ~ 10%. Such masses are quantitatively consis- 
tent with those d erived using a simple colour-stellar M/L relation 
dBell et alfc 003). and comparison of stellar and dynamical masses 
for a few z ~ 1 early-type galax ies yielded consist ent results to 
within their combined errors (see iBorch et alj|2006l for more de- 
tails). 

There are some galaxies for which the 17-band classification 
failed to find a satisfactory solution (2% of the galaxies with red- 
shift estimates); we choose to adopt in these cases a rest-frame 
colour-derived stellar mass, using rest-frame B and V absolute 
magnitudes/luminosities, and a l/-band absolute magnitude of the 
Sun of 4.82: 

log 10 M*/Mq = -0.728 + 1.305(B-y)+log 10 L V /L Q .(10) 

As with restframe photometric properties, we also present es- 
timates of stellar mass assuming that the galaxy is at the cluster 
redshift (denoted in the catalogues by the suffix '_cl' in the column 
names). Random stellar mass errors are estimated to be ~ 0.1 dex 
on a galaxy-by-galaxy basis in most cases, and systematic errors 
in the stellar masses (setting the overall mass scale and its red- 
shift evolution) were argued to be at the 0. 1 dex level for galaxies 
without ongoing or recent major starbursts; for galaxies with strong 
bursts, masses could be overestimated by < 0.5 dex. 

Finally, we note potential aperture effects on stellar masses 
and SEDs for some objects. The colours are estimated within an 
aperture but are normalized by the total light in the deep 7?-band 



9 Local comparison samples, e.g., the SDSS, typically adopt template com- 
binations with 'older' ages, potentially leading to offsets between the over- 
all mass scale of our masses and local masses at a given rest-frame colour. 
We make no attempt to resolve this issue here, and refer the interested reader 



to lBell &deJonel 120011) and lBell et all J2007n for further discussion of this 
issue. 



image alone. For small objects or particularly large objects without 
colour gradients this has no consequence. But if large size, low con- 
centration and strong colour gradients are combined, the total SED 
will deviate from the aperture SED underlying the M/L estimate. 
In a companion paper studying properties of spiral galaxies in the 
supercluster, Wolf et al. (MNRAS, accepted) have investigated this 
effect by examining the total colours across a wide parameter space 
in the sample. In most cases the aperture values are similar to the 
total ones, but they identify an issue for morphologically-classified 
spiral galaxies in the supercluster and eliminate the highest-mass 
regime with log M*/Mq > 11 from their study. 

4.4 GALEX 

The Abell 901/902 field was observed by GALEX in the far-UV 
(/, A cft ~ 1528A) and near-UV (n, A cff ~ 227lA) bands0 In- 
dividual observations (or single orbit 'visits') between the dates 
12 February 2005 and 25 Febr uary 2007 were coad ded by the 
GALEX pipeline (GR4 version iMorrissev et al.| [2007) to produce 
images with net exposure times of 57.18 ks in n (47 visits) and 
50.19 ks in / (40 visits). The GALEX field of view in both bands 
is a 0.6°radius circle, and the average centre of the visits (the 
GALEX field centre) is (a,5) J2 ooo = (9 h 56 m 20:7, -10°6'2l'/6). 
The GALEX PSF near the field centre has ~ 4.2" FWHM at / and 
~ 5.3" FWHM at n, both of which increase with distance from the 
field centre (variations in the PSF that are not a function of distance 
from the field centre are smoothed out by the distribution of roll an- 
gles of the visits). The astrometric accuracy is ~ 0.7", and > 97% 
of catalogued source positions are within 2" of their true positions. 
The phot ometric calibration is s table to 0.02 mag in n and 0.045 
mag in / dMorrissev et alj|2007h . 

Source detection and photometry is via the GALEX pipeline 
code, which employs a version of SExtractor (Bert in & Amoutsl 
modified for use with low-background images. Magnitudes 
are measured both in fixed circular apertures and in automatic Kron 
elliptical apertures, and in isophotal apertures. The 5a point-source 
sensitivities in the Abell 901/902 field are / ~ 24.7 mag (AB) and 
n ~ 25.0 mag (AB), though there are spatial variations across field, 
especially a slightly decreasing sensitivity towards the edge of the 
field. At these levels source confusion in the n band becomes an is- 
sue, and the n band fluxes of faint objects (n > 23 mag) are likely 
to be overestimated. GALEX data products include intensity, back- 
ground, and relative response (i.e., effective exposure time) maps in 
both bands as well as source catalogues in both bands and a band- 
merged source catalogue. 



4.5 2dF spectroscopy 

Spectra of cluster galaxies were obtained using the 2dF instrument 
on the AAT in March 2002 and March 2003. A total of 86 galax- 
ies were observed using the 1200B grating (spanning the observed 
wavelength range 4000-5100 A) in a single fibre configuration dur- 
ing the 2002 run. Three fibre configurations using the lower resolu- 
tion 600V grating (spanning 3800-5800 A) were observed during 
the 2003 run: fibres were placed on 368 objects, with 47 repeated 
from 2002. The primary selection function assigned higher prior- 
ity to those galaxies selected by photometric redshift to be within 



10 Unlike all other datasets detailed here, the GALEX observations were 
not led by members of the STAGES team. We list the publicly archived data 
products here for completeness. 
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the supercluster redshift slice and having R < 20, with additional 
fibres being allocated to secondary targets (including fainter galax- 
ies and a small number of white dwarfs and QSOs) when available. 
Data reduction was performed with the standard 2dfdr (v2.3) 
pipeline package. 

In total, spectra were obtained for 407 unique objects. Red- 
shifts were determined by two independent means: firstly by man- 
ual line profile fitting of the Ca H and K features in absorp- 
tion and secondly by cross-correlation with template spectra us- 
ing the XCSAO task within IRAF (Kurtz & Mink 1998). Compar- 
ison of the two measurements showed no cause for concern, with 
A z = 0.00149 ± 0.00006. After eliminating non-galaxy and poor 
quality spectra, we have redshifts for 353 galaxies in total. 

The 2dF spectroscopic data have previously been used to 
quantify the reliability of the COMBO-17 redshifts in W04 (see 
also j{3]>, to verify cluster me mbership for the matched X-ray point 
sources ( Gilmour et al. 2007)), and to create composite spectra for 
three photometric classes of cluster galaxies in WGM05. A dynam- 
ical analysis of the the clusters using the 2dF redshifts will be pre- 
sented in Gray et al. (in prep.). 



4.6 XMM-Newton 

X-ray data for the A901/2 region is desirous both to detect point- 
source emission from cluster members (star-formation or AGN) 
and the extended intracluster medium (ICM). A 90 ks XMM im- 
age of the A901/2 field was taken on May 6/7 2003 using the three 
EPIC cameras (MOS1, MOS2 and PN) and a thin filter, under pro- 
gram 14817 (PI: Gray). The level 1 data were taken from the sup- 
plied pipeline products, and reduced with SAS v5.4 and the calibra- 
tion files available in May 2003. Final exposure times were ~ 67 
ks for MOS and ~61 ks for PN following the removal of time inter- 
vals suffering from soft proton flares. Four energy bands were used: 
0.5-2 keV (soft band), 2-4.5 keV (medium band), 4.5-7.5 keV (hard 
band) and 0.5-7.5 keV (full band). 

The creation of the point-source catalogue usin g wavelet de- 
tectio n methods is described in detail elsewhere l lGilmour et al.l 
120071) . A total of 139 significant sources were found. The pres- 
ence of an X-ray luminous Type-I AGN near the centre of A901a 
(see Appendix IA3t complicated the detection of the underlying 
extended cluster emission. A maximum-likelihood technique was 
used to match this catalogue to COMB O-17 resulting i n 66 s ecure 
counterparts with photometric redshifts. iGilmour et al.l j2007l) used 
these data to examine the local environments of the cluster AGN 
and their host properties. 

To isolate the remaining extended emission coming from the 
clusters, a separate conservative point-source catalogue was con- 
structed. Care was taken to remove both the cosmic background and 
spatial variations in the non-cosmic background. The background 
subtracted images were weighted by appropriate energy conversion 
factors to create flux images for each detector. These flux images 
were masked and summed together to create merged background- 
subtracted images in each band. 

Point source regions were removed and replaced with the local 
background value selected randomly from a source free area within 
10 pixels (or 20 pixels if there were not enough background pixels 
within the smaller radius). Smoothed images were created in each 
band using a Gaussian kernel of radius 4 pixels. Maps of the ex- 
tended emission and an examination of the global X-ray properties 
of the clusters will be presented in Gray et al. (in prep.). 



4.7 GMRT 

The A901/2 field was observed on 2007 March 25th and 26th 
March with the Giant Metrewave Radio Telescope (GMRT, see 
lAnanth akrishnan 2005 for further details). The field was centred 
at (a,(5)j2ooo = (09 h 56 m 17 s , -10°01'28") and observed at 610 
and 1280 MHz on respective nights. The GMRT is an interferome- 
ter, consisting of thirty antennas, each 45 m in diameter. The bright 
sources 3C147 and 3C286 were observed at the start and end of 
each observing session, in order to set the flux density scale. During 
the observations a nearby compact source 0943—083 was observed 
for about 4 minutes at roughly 30 minutes intervals, to monitor and 
correct any antenna-based amplitude and phase variations. 

The total integration time on the field was ~6.5 hours at each 
frequency. The observations covered two 16 MHz sidebands, posi- 
tioned above and below the central frequency. Each sideband was 
observed with 128 narrow channels, in order to allow narrow band 
interference to be identified and efficiently removed. The observed 
visibility data were edited and calibrated using standard tasks with 
the AIPS package, and then groups of ten adjacent channels were 
averaged together, with some end channels discarded. This reduced 
the volume of the visibility data, whilst retaining enou gh channels 
so tha t chromatic aberration is not a problem (e.g., see lGarn et al.l 
2007, for further details of GMRT analysis). Given the relatively 
large field of view of the GMRT compared with its resolution, 
imaging in AIPS requires several 'facets' to be imaged simulta- 
neously, and then be combined. Preliminary imaging results, af- 
ter several iterations of self-calibration, have produced images with 
resolutions of about 5"and 2'.'5 at 610 and 1280-MHz respectively, 
with r.m.s. noises of approximately 25 and 20 /iJy beam" 1 in the 
centre of the fields, before correction for the primary beam of the 
GMRT. The primary beam - i.e. the decreasing sensitivity away 
from the field centres due to sensitivity of individual 45-m anten- 
nas - is approximately Gaussian, with a half-power beam width 
(HPBW) of approximately 44' and 26' at 610 and 1280-MHz re- 
spectively. These images are among the deepest images made at 
these frequencies with the GMRT. Further analysis and the source 
catalogue will be presented in Green et al.(in prep.). 

4.8 Simulations and mock galaxy catalogues 

In order to facilitate the interpretation of the observational re- 
sults and to study the physical processes of galaxy evolution, N- 
body, hydrodynamic, and semi-analytic simulations that closely 
mimic the A901/2 system are being produced (van Kampen et 
al., in prep.). We c onstrain initial conditions using the method of 
Hoffman & Ribak] Jl99lh to take into account the gross properties 
of A901a, A901b, A902, the SW group, and the neighbouring clus- 
ters A868 and A907 (outside the observed field). The simulations 
produce a range of mock large-scale structures to test three ba- 
sic formation scenarios: a 'stationary' case, where A901(a,b) and 
A902 will not merge within a Hubble time, and a pre- as well as a 
post-merger scenario. When the likelihood of each scenario is un- 
derstood, one can further test the models for the detailed physical 
processes known to be operating on galaxies in and around such 
clusters. 



5 SUMMARY AND DATA ACCESS 

We have presented the multiwavelength data available for the 
A901/2 supercluster field as part of the STAGES survey: high- 
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resolution HST imaging over a wide area, extensive photometric 
redshifts from COMBO-17, and further multiwavelength observa- 
tions from X-ray to radio. These data have already been used to 
create a high r esolution mass map of the system using weak gravita- 
tional lensing dHevmans et al .120081) . Further work by the STAGES 
team to study galaxy evolution and environment is ongoing and in- 
cludes the following: 

• Gallazzi et al. I J2008h explore the amount of obscured star- 
formation as a function of environment in the A901/2 superclus- 
ter and associated field sample by combining the UV/optical SED 
from COMBO-17 with the Spitzer 24/im photometry in galaxies 
with M* > 10 10 Mq. Results indicate that while there is an over- 
all suppression in the fraction of star-forming galaxies with density, 
the small amount of star formation surviving the cluster environ- 
ment is to a large extent obscured. 

• Wolf et al. (MNRAS, accepted) investigate the properties of 
optically passive spiral and dusty red galaxies in the superclus- 
ter and find that the two samples are largely equivalent. These 
galaxies form stars at a substantial rate that is only a factor of 
four times lower than blue spirals at fixed mass, but their star for- 
mation is more obscured and has weak optical signatures. They 
constitute over half of the star forming galaxies at masses above 
log M» /Mq — 10 and are thus a vital ingredient for understand- 
ing the overall picture of star-formation quenching in cluster envi- 
ronments. 

• Marinova et al. (ApJ, submitted) identify and characterize bars 
in bright (My ^ —18) cluster galaxies through ellipse-fitting. 
The selection of moderately inclined disk galaxies via three com- 
monly used methods, visual classification, colour, and Sersic cuts, 
shows that the latter two methods fail to pick up many red, bulge- 
dominated disk galaxies in the clusters. However, all three meth- 
ods of disk selection yields a similar global optical bar fractions 
(/bar-opt ~ 0.3), averaged over all galaxy types. When host 
galaxy properties are considered, the optical bar fraction is found 
to be a strong function of both the luminosity and morphological 
property (bulge-to-disk ratio) of the host galaxy, similar to trends 
recently reported in field galaxies. Furthermore, results indicate 
that the global optical bar fraction for bright galaxies is not a strong 
function of local environment. 

• Heiderman et al. (in prep.) identify interacting galaxies in the 
supercluster using quantitative analysis and visual classifications. 
Their findings include that 4.9 ± 1.3% of bright (M v sC -18), in- 
termediate mass (M* ^ 1 x 10 9 Mq) galaxies are interacting. The 
interacting galaxies are found to lie outside the cluster cores and to 
be concentrated in the region between the cores and virial radii of 
the clusters. Explanations for the observed distribution include the 
large galaxy velocity dispersion in the cluster cores and the possi- 
bility that the outer parts of the clusters are accreting groups, which 
are predicted to show a high probability for mergers and strong in- 
teractions. The average star formation rate is enhanced only by a 
modest factor in interacting galaxies compared to non -interacting 
galaxi es, similar to conclusions reported in the field by jjogee et al.l 
( 2008). Interacting galaxies only contribute ~ 20% of the total SFR 
density in the A901/902 clusters. 

• Boehm et al. (in prep.) are utilizing the stability of the PSF on 
the STAGES images for a morphological comparison between the 
hosts of 20 type-1 AGN and 200 inactive galaxies at an average 
redshift (z) ~ 0.7. This analysis includes extensive simulations 
of the impact of a bright optical nucleus on quantitative galaxy 
morphologies in terms of the CAS indices and Gini/M20 space. 
We find that the majority of the hosts cover parameters typical for 



disk+bulge systems and mildly disturbed galaxies, while evidence 
for strong gravitational interactions is scarce. 

• Bacon et al. (in prep.) are examining the higher order lens- 
ing propert ies of the STA GES data. They construct a shapelets 
catalogue teefregiedl2003h for the STA GES galaxies; this is then 
used to estimate the gravitational flexion dBacon et al J 20061) at each 
galaxy position. Galaxy-galaxy flexion is measured, leading to esti- 
mates of concentration and mass for STAGES galaxies; constraints 
on cosmic flexion are also found, showing very good containment 
of systematic effects. The ability of flexion to improve convergence 
maps is also discussed. 

• Robaina et al. (in prep.) make use of a combined GEMS and 
STAGES sample of 0.4 < z < 0.8 galaxies to find that interacting 
and merging close pairs of massive galaxies (> 1O 1O M0) show a 
modest enhancement of their star formation rate; in particular, less 
than 15% of star formation at 0.4 < z < 0.8 is triggered by major 
interactions and mergers. 

• Barden et al. (in prep.) are exploring both the GEMS and 
STAGES data sets to investigate the evolution of structural param- 
eters of disc galaxies as a function of luminosity and stellar mass 
over a wide range of environments and morphologies. In the pro- 
cess, GALAPAGOS will be extended to perform bulge/disc decom- 
position. 

• Mcintosh et al. (in prep.) are using both quantitative and qual- 
itative morphologies to explore the morphological mix of red se- 
quence galaxies as a function of stellar mass over the last seven 
billion years from the combined STAGES + GEMS sample. 

It is our intention that the data products described here should 
be publicly available for use by the wider community for those 
interested in the supercluster itself or for data-mining the entire 
survey volume. To that end, the reduced HST images (both tiles 
and individual galaxy postage stamps) are available for download 
at the Multimission Archive at Space Telescopf3 (MAST). Fur- 
thermore, the complete STAGES catalogue described in this pa- 
per is available from the STAGES websiteQ including all HST- 
derived parameters; GALFIT profile fitting results; COMBO-17 
photometry, SEDs and photometric redshifts; and stellar masses 
and star-formation rates. The multiwavelength data available there 
includes the Spitzer/MI PS 24[i.m images and catalogue; the X-ray 
point source catalogue ( Gilmour et al. 2007) and the gravitational 
lensing mass maps ( IHevmans et alj|2008l) . GALEX data and cata- 
logues are available via MAST. The X-ray maps, 2dF spectra and 
radio catalogue and mocks will be also be placed on the website 
with the publication of their associated papers, or may be made 
available upon request. Table [8] contains a summary of the avail- 
able data products. 
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Table 8. Description of all available A901/902 data products. 
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immediate 
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STAGES master catalogue: SExtractor, GALFIT, COMBO-17, stellar masses, SFRs 
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this paper 


LUMdU-1 / abUj ana completeness tables 
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this paper 


GALM1 pronle fitting completeness from simulations 
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this paper 


Spitzer 24 |am imaging and catalogue 


immediate 


this paper 


HST-derived weak lensing mass map 


immediate 


Heymans et al 2008 


XMM point source catalogue 


immediate 


Gilmour et al 2007 


GALEX imaging and catalogues (from the GALEX archive) 


immediate 


this paper 


X-ray imaging 


on request 


Gray et al. (in prep.) 


2dF spectroscopy 
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GMRT catalogue 


TBC 


Green et al. (in prep.) 
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Figure Al. The COMBO-17 SEDs of the merging system 44635 (top) and 
45154 {bottom). The latter case is a dust-reddened fit by eye to z = 0.08, 
the likely redshift of the system. 
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APPENDIX A: NOTES ON INDIVIDUAL OBJECTS 

Here, we collect some details on ten noteworthy individual ob- 
jects that have either extreme properties or are intrinsically rare and 
found only by chance in a field of this size. They are drawn from 
the COMBO-17 sample and are identified here via their COMBO- 
17 object numbers. 

Al The brightest near-infrared source: a Mira variable 

The object with the COMBO-17 number 35250 is classified as a 
very red star of spectral type M8 III in the 13th Genera l Cat- 
alogue of MK Spectral Classifications dBuscombel [l998). It is 
also known as the IRAS point source 09540 — 0946, located at 
(a,5)j20oo = (09 fa 56 m 32^4 -10 o 0l'l5"), an d it is a ROSAT 
All-Sky Survey Bright source jVogesetal.ll 19991) . It has a very red 
SED with {B, R, J, K) « (16, 13.3, 7.25, 5.75) and is the bright- 
est object in the field at A > 1 \im. However, it has a large variabil- 
ity amplitude and was identified as a long-period pulsating Mira 




Figure A2. Left panel: The 20 ks COMBO-17 i?-band image of the merg- 
ing system that is with ~ 50 mjy the brightest extragalactic 2Ap source 
in the field (objects 44635 and 45154, size of image 1' X 1', N is up, E is 
left) and missing from the matched catalogue. Right panel: The same image 
in hard cuts reveals a tidal arm with l/5000th of the surface brightness of 
the central disks. This arm is too faint to be visible in the STAGES/HST 
images. 




Figure A4. The COMBO-17 SED of object 12716, the central dominating 
(cD) galaxy of CB I, the cluster at z 0.47 in the background of A902. 

star in a search for high-redshift QSOs (Kirkp atrick et alJI 1997h . 
The area around this object had to be excluded from the Spitzer 
IRAC imaging due to its high brightness. 

A2 The brightest far-infrared galaxy: a merger 

The brightest 24[im galaxy is a system of two merging disk galax- 
ies with a total magnitude of R ~ 16. The Northern system (44635) 
has a very blue SED (U — V) T est = 0.14 and implies a very strong 
Ha line given its elevated i?-band flux (see Fig. lAlt . The Southern 
system (45154) has an extremely red SED (U - V) rcs t = 1.97 
and implies strong dust-reddening. Their reshifts are estimated as 
-Zphot — 0.084 and 0.053, but the blue SED is better constrained by 
emission lines. Assuming z — 0.08 for both objects, the projected 
separation between their two nuclei of 3'.' 5 translates into 5 kpc. 

The 20 ks 7?-band image of COMBO-17 shows tidal features 
with very low surface brightness (Fig. IA2b . The arm that reaches 
once around the entire galaxy has 5000 x lower surface brightness 




Figure A5. The COMBO-17 SED of STAGES I (object 59586), a dwarf 
irregular at z phot ss 0.04 (but likely z < 0.01). 
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Figure A6. The COMBO-17 SED of the only photometrically-classified 
'strange' object in the dataset: galaxy 54511 isatz 0.3 and has extremely 
strong emission lines (OIII+H/3 with EW ~ 150 nm). 




Figure A7. The COMBO-17 SED of object 474, the bluest white dwarf in 
the field. The lack of H/3 absorption (see 485 filter) makes it a DB white 
dwarf. The best-fitting temperature is ~ 30, 000 K. 



than the main disks of the two merging galaxies. The system is 
also a strong radio source (NVSS J095643-095544) and was seen 
by IRAS. In our Spitzer MIPS images it shows ~ 50 mjy of flux, 
but such bright FIR measurements are missing from our matched 
catalogue due to matching difficulties. Preliminary analysis of the 
GMRT data reveals a strong radio detection at both 1280 MHz 
and 610 MHz with total flux S(1280)=5.63 ± 0.05 mjy/bm and 




Figure A8. The COMBO-17 SED of object 33783, the faintest white dwarf 
in the field. The strong H/3 absorption line in the 485 filter allows its classi- 
fication even at this faint level (R = 23.4, T cB ss 11, 000 K). 



S(610)=13.68±0.05 mjy. The radio souce is partially resolved with 
a deconvolved size of 4.5" x 1.9" at 1280 MHz and 4.1" x 1.9" 
at 610 MHz, and a position angle of 40° at both frequencies. 

A3 The brightest X-ray source: a type-I AGN in A901a 

Object 41435 is a massive red -sequence elliptical w ith excess blue 
light in its SED (see Fig. 7 of lGilmour et al.ll2007h that has biased 
the redshift estimation. While it has z p hot ~ 0.33, it is almost 
certainly a cluster member and a z = 0.16 template fitted by hand 
works well and leaves over some room for AGN light. It is the 
brightest X-ray source in the STAGES field observed by XMM and 
a point source with a luminosity (assuming z — 0.16) of Lx = 
1.55 x 10 44 erg/s. It is also the brightest radio source at 1280 MHz 
and is unresolved with total flux 46.33±0.01 mjy/bm. At 610 MHz 
it is partially resolved with an integrated flux density of 171.3±0.1 
mjy. 



A4 An SO galaxy with a full Einstein ring 



Object 14049 (Fig. lA3l left) is an SO galaxy displaying a full optical 
Einstein ring. It has z p hot = 0.23, but 2dF spectroscopy confirms 
it is a cluster member with z spcc — 0.168, implying that the SED 
is contaminated by light from the lensed galaxy. Subsequent tar- 
getted spectroscopy revealed a source redshift z 3 = 1.5 (Aragon- 
Salamanca et al, in prep.). 

A5 A galaxy cluster in projection behind A902: CBI 

Examination of the redshift distribution along the line-of-sight to 
the A902 cluster revealed the presence of a massive background 
cluster at z ~ 0.47, subseq uently designated CBI (Fig. [A3] centre). 
A 3D lensing approach dTavlor et al.1 2004) was used to constrain 
the m asses of the two c lusters beyond the 2D mass reconstruc- 
tion of iGrav et al] J2002I) . Object 12716 (R = 19.1) is the central 
cD galaxy of CBI and is detected as an unresolved object in the 
preliminary analysis of the GMRT data with S int (1280MHz) = 
2.07 ± 0.02mJy/bm and SU(610MHz) = 6.70 ± 0.04 mjy. Its 
brighter and bluer close neighbour (R — 17.8) is an actual member 
of A902 (see Fig.|A3] centre). 



A6 The dwarf irregular galaxy STAGES I 

The object with the COMBO-17 number 59586 is a nearby dwarf 
irregular galaxy (see Fig. I A3 1 right and Fig. IA5t estimated at 
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-Zphot = 0.044 ± 0.026 (consistent with z = at 1.6cr). At the 
estimated redshift it would have My ~ —16.7 and log M*/Mq « 
8.7; however given the brightness of the resolved point sources it 
is most likely at z < 0.01. It has a Sersic index of n — 0.55 and 
shows clear signs of irregularity besides a blue colour. 

A7 The galaxy with the strongest emission lines 

The COMBO- 17 catalogue contains only one object classified as 
'strange' as a result of having a Xrcd > 30 for its best template fit, 
while having good flags: object 54511 is a galaxy with extremely 
strong emission lines and R ~ 22.5. The emission-line flux in the 
_R-band and the 646-band both suggest EW « 150 nm, which 
would need to be the combined H/3 and OIII lines. A line in the 
filter 485/30 shows EW w 14 nm and is possibly OIL The redshift 
of the object appears to be constrained to 0.27 < zn ncs < 0.32 by 
a third line signal in the filter 855 (Ha; see Fig. lA6b . 

A8 The bluest white dwarf: U — B < —1 

Object 474 is the bluest white dwarf with a satisfying fit to our 
DA template library, although the SED (see Fig. |A7| > shows clearly 
no H/3 absorption line, rendering this object a DB. The best-fitting 
temperature is ~ 30, 000 K. 

A9 The faintest white dwarf we could identify 

Object number 33783 is the faintest white dwarf our classification 
can identify with U = 23.3 and R — 23 .4. At this magnitude level, 
the WD selection is already highly incomplete, but the strong H/3 
absorption still constrains the template fit (see Fig. lA8b . 



APPENDIX B: STAGES MASTER CATALOGUE 

The tables in this section contain information relating to the 
publically-available STAGES master catalogue. Table IB ll lists and 
defines the column names containing STAGES and COMBO-17 
data and derived stellar masses and star-formation rates. Table IB2l 
details the three sample flags in the catalogue and describes how 
they are to be used to select relevant populations from the overlap 
between the HST, COMBO-17, and Spitzer datasets. 

This paper has been typeset from a TpX/ ETpX file prepared by the 
author. 



STAGES: the Space Telescope A901/2 Galaxy Evolution Survey 27 



Table Bl. Column entries in the published FITS catalogue, their headers and meanings. Some restframe luminosities are extrapolated in some redshift ranges. 
We give the redshift intervals, where no extrapolation errors are expected. 



STAGES information 


sLnumber 


object number 


st_x_image 


x-position from SExtr in [pix] on tile 


st_y .image 


y-position from SExtr in [pix] on tile 


st_cxx_image 


ellipse parameter from SExtr in [pix] 


st_cyy .image 


ellipse parameter from SExtr in [pix] 


sLcxy .image 


ellipse parameter from SExtr in [pix] 


st_theta_image 


pos. angle from SExtr in [deg] in image 




coordinates (measured from right to up) 


st_theta_world 


pos. angle in [deg] in world coordinates 


st_ellipticity 


ellipticity from SExtr 


st_kron_radius 


Kron radius in units of [st_a_image] 


st_a_image 


semi-major half-axis from SExtr in [pix] 


st_b_image 


semi-minor half-axis from SExtr in [pix] 


st_alpha_J2000 


right ascension from SExtr in [deg] 


st_delta_J2000 


declination from SExtr in [deg] 


st_ background 


background value from SExtr in [counts] 


st_flux_best 


"best" flux from SExtr in [counts] 


st_fluxeiT_best 


eiTor of st_flux_best 


st_mag_best 


"best" magnitude from SExtr in [AB mag] 


st_magerr_best 


error of st_mag_best 


st_flux_radius 


half-light radius from SExtr in [pix] 


st_isoarea_image 


isophotal area from SExtr in [pix 2 ] 


st_fwhmjmage 


FWHM from SExtr in [pix] 


sLflags 


SExtr quality flags 


st_class_star 


SExtr stellarity estimator 


st_org_image 


postage stamp image file name 


st_file_galfit 


GALFIT output filename containing fit data 


st.X.galfit 


x-position on postage stamp in [pix] 


sLXerr.galfit 


error of st_X_galfit 


st_Y_galfit 


y-position from GALFIT in [pix] 


st_Yerr_galfit 


eiTor of st_Y_galfit 


st.MAG.galfit 


total magnitude from GALFIT in [AB mag] 


st_MAGerr_galfit eixor of st_MAG_galfit 


st_RE_galfit 


half-light radius from GALFIT in [pix] 


st_REerr_galfit 


eiTor of st_RE_galfit 


stJM_galfit 


Sersic index from GALFIT 


st_Nerr_galfit 


error of st_N_galfit 


st_Q_galfit 


major-to-minor axis ratio from GALFIT 


st_Qerr_galfit 


eiTor of st_Q_galfit 


st.PA.galfit 


pos. angle in [deg] measured from up to left 


st_PAerr_galfit 


error of st_PA_galfit 


st_sky_galfit 


sky value from GALAPAGOS 


st.tile 


tile number in STAGES mosaic 


COMBO-17 general information 


COMBO-nr 


COMBO-17 A901/2 field object number 


ra 


right ascension (J2000) 


dec 


declination (J2000) 


xpix 


x-position on COMBO-17 f?-frame in pixels 


ypix 


y-position on COMBO-17 f2-frame in pixels 


Rmag 


total i?-band magnitude 


e_Rmag 


1-cr error of total i?-band mag 


ap_Rmag 


aperture f?-band magnitude in run E 


apd_Rmag 


difference total to aperture (point source ~ 0) 


Various flags for sample selection 


photjlag 


COMBO-17 photometry flags (see Sect. 3.5) 


combo.flag 


COMBO-17 sample flag (see TablelB2l 


stages_flag 


STAGES sample flag (see TablelB2l 


mips_flag 


MIPS sample flag (see TablelB2l 



COMBO-17 classification results 



chi2red 

chi2reds 

chi2redg 

chi2redq 

chi2redw 

chi2redg_cl 

mc.class 

mc_z 

e_mc_z 

mc_z2 

e_mc_z2 

mc_z_ml 

mc_Ebmv 

e_mc_Ebmv 

mc_Ebmv_ml 

mc_age 

e_mc_age 

mc_age_ml 

mc_z_cl 

mc_Ebmv_cl 

e_mc_Ebmv_cl 

mc_age_cl 

e_mc_age_cl 



X 2 /Nf of best-fitting template 

X 2 /Nf of best-fitting star template 

X 2 /Nf of best-fitting galaxy template 

X 2 /-/V/ of best-fitting QSO template 

X 2 /Nf of best-fitting WD template 

X 2 /Nf of best-fitting galaxy template at z = 0.167 

multi-colour class (see Table[6} 

mean redshift in distribution p(z) 

standard deviation (1-cr) in distribution p(z) 

alternative redshift if p(z) bimodal 

standard deviation ( 1 -cr) at alternative redshift 

peak redshift in distribution p(z) 

mean E(B — V) in distribution p(z) 

standard deviation (1-cr) in distribution p(E(B-V)) 

peak value in distribution p(E(B — V)) 

mean template age index 

standard deviation (1-cr) of template age index 

peak in template age index distribution 

redshift assuming cluster membership 

mean E(B — V) assuming cluster membership 

standard deviation in p(E(B — V)) if cluster member 

mean age index assuming cluster membership 

standard deviation in age index if cluster member 



total galaxy restframe luminosities 



S280Mag 


A/ab s ,gai in 280/40 (z w [0.25, 1.3]) 


e_S280Mag 


1-cr error of M abSjgal in 280/40 


UjMag 


M a b s , g ai in Johnson U (ok at all z) 


e.UjMag 


1-cr error of Af a bs gal i n Johnson U 


BjMag 


A^abs,gal in Johnson B (z »* [0.0, 1.1]) 


e_BjMag 


1-cr eiTor of M a bs,gal i n Johnson B 


VjMag 


Af a bs,gal in Johnson V (z a [0.0, 0.7]) 


e.VjMag 


1-cr error of Af a bs gal in Johnson V 


usMag 


A^abs gal in SDSS u (ok at all z) 


e.usMag 


1-cr eiTor of Af a b 8 g a i in SDSS u 


gsMag 


Afabs.galinSDSSsO^ [0.0,1.0]) 


e.gsMag 


1-cr error of M a b s ,gal in SDSS g 


rsMag 


Afabs.gal in SDSS r (z [0.0, 0.5]) 


e_rsMag 


1-cr eiTor of M a bs,gal in SDSS r 


restframe luminosities at cluster distance 


S280Mag.cl 


Afabs, gal in 280/40 (if cluster member) 


e.S280Mag.cl 


1-cr error of M a b s ,gal in 280/40 


UjMag.cl 


Afabs, gal in Johnson U (if cluster member) 


e.UjMag.cl 


1-cr error of Af a bs gal in Johnson U 


BjMag.cl 


Afabs, gal in Johnson B (if cluster member) 


e_BjMag_cl 


1-cr error of Af a bs gal in Johnson B 


VjMag.cl 


Afabs, gal in Johnson V (if cluster member) 


e.VjMag.cl 


1-cr error of Af a bs gal in Johnson V 


usMag.cl 


Afabs, gal in SDSS u (if cluster member) 


e_usMag_cl 


1-cr error of Af a b s ,gal in SDSS u 


gsMag.cl 


Afabs, gal in SDSS g (if cluster member) 


e_gsMag_cl 


1-cr error of M a b s ,gal in SDSS g 


rsMag.cl 


M a bs,gal in SDSS r (if cluster member) 


e_rsMag_cl 


1-cr error of Af a bs,gal in SDSS r 


QSO restframe luminosities 


S145Mag 


Afabs.QSO in 145/10 (z [1.4, 5.2]) 


e.S145Mag 


1-cr error of Af ab s,QSO in 145/10 
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Table Bl - continued 



observed seeing-adaptive aperture fluxes 



W420f 


photon flux in filter 420 


e.W420f 


1-<t photon flux error in 420 


W462f 


photon flux in filter 462 


e.W462f 


i-cr photon flux error in 462 


W485f 


photon flux in filter 485 


e.W485f 


1-cr photon flux error in 485 


W518f 


photon flux in filter 518 


e.W518f 


i-cr photon flux error in 518 


W571f 


photon flux in filter 571 


e.W571f 


1-<t photon flux error in 571 


W604f 


photon flux in filter 604 


e.W604f 


i-cr photon flux error in 604 


W646f 


photon flux in filter 646 


e.W646f 


i-cr photon flux error in 646 


W696f 


photon flux in filter 696 


e.W696f 


i-cr photon flux error in 696 


W753f 


photon flux in filter 753 


e.W753f 


i-cr photon flux error in 753 


W815f 


photon flux in filter 815 


e.W815f 


i-<j photon flux error in 815 


W856f 


photon flux in filter 856 


e.W856f 


i-cr photon flux error in 856 


W914f 


photon flux in filter 914 


e.W914f 


i-cr photon flux error in 914 


Uf 


photon flux in filter U 


e.Uf 


i-cr photon flux error in U 


Bf_A 


photon flux in filter B in run A 


e_Bf.A 


i-cr photon flux eiTor in BIA 


Bf.G 


photon flux in filter B in run G 


e_Bf_G 


i-cr photon flux error in BIG 


Vf 


photon flux in filter V 


e.Vf 


i-<j photon flux error in V 


Rf 


photon flux in filter R 


e_Rf 


i-<j photon flux error in R 


If 


photon flux in filter I 


ejf 


i-cr photon flux error in I 


observed aperture Asinh Vega magnitudes 


W420magA 


magnitude in filter 420 


e.W420magA 


i-cr magnitude error in 420 


W462magA 


magnitude in filter 462 


e_W462magA 


i-cr magnitude error in 462 


W485magA 


magnitude in filter 485 


e_W485magA 


i-cr magnitude error in 485 


W518magA 


magnitude in filter 5 1 8 


e.W518magA 


i-cr magnitude error in 518 


W571magA 


magnitude in filter 571 


e.W571magA 


i-cr magnitude error in 571 


W604magA 


magnitude in filter 604 


e_W6()4magA 


i-cr magnitude error in 604 


W646magA 


magnitude in filter 646 


e_W646magA 


i-cr magnitude error in 646 


W696magA 


magnitude in filter 696 


e_W696magA 


i-cr magnitude error in 696 


W753magA 


magnitude in filter 753 


e.W753magA 


i-cr magnitude error in 753 


W815magA 


magnitude in filter 815 


e.W815magA 


i-cr magnitude error in 815 


W856magA 


magnitude in filter 856 


e.W856magA 


i-cr magnitude error in 856 


W914magA 


magnitude in filter 914 


e_W914magA 


i-cr magnitude error in 914 



observed aperture Asinh Vega magnitudes (cont.) 



UmagA 


magnitude in filter U 


e.UmagA 


I-cr magnitude error in U 


BmagA_A 


magnitude in filter B in run A 


e_BmagA_A 


i-a magnitude error in BIA 


BmagA.G 


magnitude in filter B in run G 


e_BmagA_G 


i-a magnitude error in BIG 


VmagA 


magnitude in filter V 


e.VmagA 


l-cr magnitude error in V 


RmagA 


magnitude in filter _R 


e_RmagA 


l-<7 magnitude error in R 


ImagA 


magnitude in filter / 


eJmagA 


1-cr magnitude error in / 


stellar masses and star formation rates 


logmass 


log 10 of stellar mass 


logmass_cl 


log 10 of stellar mass if cluster member 


flux24 


MIPS 24// flux in microjy 


tir 


IR luminosity in Lq 


tuv 


UV luminosity in Lq 


tir.cl 


IR luminosity in Lq if cluster member 


tuv.cl 


UV luminosity in Lq if cluster member 


sfr.det 


SFR from UV + IR if IR detected 


sfrjo 


SFR lower limit from UV alone 




(if IR non-detected) 


sfr_hi 


SFR upper limit (if IR non-detected) 


sfr_det_cl 


SFR if IR detected (if cluster member) 


sfr_lo_cl 


SFR lower limit from UV alone 




(if no-IR, if cluster member) 


sfr_hi_cl 


SFR upper limit (if no-IR, if cluster member) 


sed.type 


l=old red, 2=dusty red, 3=blue cloud 


sed_type_cl 


l=old red, 2=dusty red, 3=blue cloud (if cluster member) 
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Table B2. Sample flags in the public FITS catalogue and their meaning. Note that due to a m anual reinspection of COMBO-17 photometric quality flags for 
this work, the 'WGM05' sample contains 9 fewer objects than the actual published sample of lWolf et alj feOOSl) . However, we retain the name for simplicity. 
As an example, to select objects that are defined by COMBO-17 photometry as galaxies and also have extended morphologies on the HST imaging, one would 
require that combo.flag 3 and stages_fiag ^ 3. 



i lag 


Value 


Definition 


N 


STAGES .FLAG 





not in STAGES footprint (only in COMBO-17) 


6577 




1 


in STAGES footprint, but not detected by STAGES (only in COMBO-17) 


6497 




2 


detected by S lAGhS, but not HS I extended source 


5061 




3 


HST extended source but GALFIT run into constraint 


16123 




4 


HST extended source but GALFIT successful 


54621 


COMBOJTAG 





not in COMBO-17 footprint (only in STAGES) 


1271 




1 


in LUMbU-1 / tootpnnt, but not detected by LUMdU-1 / (only in ^ 1AUEIS) 


23833 




2 


detected by COMBO-17, but neither galaxy, nor cluster, nor WGM05 


48860 




3 


galaxy but neither cluster, nor WGM05 


12625 




4 


cluster galaxy, but not WGM05 


1504 




5 


cluster galaxy in WGM05 


786 


MIPSJTAG 





detected only by STAGES 


25104 




1 


detected by COMBO-17, but outside MIPS footprint 


11858 




2 


detected by COMBO-17 and inside MIPS footprint, but not detected by MIPS 


48885 




3 


detected by COMBO-17 and detected by MIPS 


3032 
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